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Summary 


A list is given of the spectral components (Types II, III and V bursts, 
and storms with continuum) of all major solar radio events recorded with the 
Dapto radio spectrograph (25-210 Mc/s) between August 1952 and September 
1960. Where possible the associated flare and ionospheric data are appended. 

The list shows that during the period in question the flares associated with 
major disturbances were twice as numerous in the northern hemisphere of the 
Sun than in the southern. The list also suggests that the probability of a 
Type II burst being followed by a continuum storm is much less during the 
early part of the sunspot cycle than in the later part. 





In its most complete form the major solar-flare event is seen on metre- 
wavelength spectral records to consist of several components (1). The first phase, 
which usually begins between the start and the maximum of the Ha flare and 
lasts for a few minutes (see Fig. 1), consists of a group of Type III (fast-drift) 
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Fic. 1.—Idealized sketch of spectral record showing the spectral components of a complete major 
radio event. 


bursts sometimes followed by a brief continuum burst of ‘Type V. After a lull 
of a further few minutes a second phase begins with the onset of a slow-drift 
Type IL burst which continues for the order of ten minutes. There may then 
follow a storm of long duration (e.g. hours) consisting of a high level of continuum 
radiation, sometimes with superposed Type I bursts. In general, one or more 
of the spectral components may appear without the remainder. Thus by far 
the commonest type of flare-associated radio disturbance is that containing a 
first-phase event only. The major events marked by Type II bursts or starts 
of continuum storms are much rarer, and are of special significance owing to 
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The spectral components of 150 major solar radio events 


No. 2, 1961 
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their relations with other solar geomagnetic and cosmic-ray effects. It is the 
purpose of the present communication to list the spectral components of all 
such major events recorded with the Dapto (N.S.W.) radio spectrograph 
between 1952 and 1960, and to give also data on their related optical flares and 
ionospheric effects. This listing is intended primarily for the use of workers 
in allied fields, but a few comments on statistical features of the table are 
appended. 

The following remarks explain the use of the table: 

Dates and times are in U.T. ‘Total times of observation are not given but 
may be found in the Quarterly Bulletin of Solar Activity. 

The radio events refer to the frequency range 40-210 Mc/s up to 1958 
November, and 25-210 Mc/s thenceforth. They are described with the aid 
of the notation used in the Ouarterly Bulletin, viz. g and G signify groups of 
Type II1 bursts containing fewer than and more than ro bursts respectively, 
and V signifies that the Type III event includes a Type V continuum event. 
A subjective indication of intensity or importance is denoted by three different 
faces of type: roman, italics and bold in ascending order. 

Entries in the column ‘‘ Storm with continuum ”’ are often referred to as 
Type IV events; however they should not necessarily be taken to indicate the 
presence of an extended, outward-moving source implied by Boischot’s original 
definition of Type IV, since the recognition of this phenomenon requires 
additional (directional) information; nor need they necessarily be related to 
events in the centimetre-wave spectrum, which have also been called Type IV. 
It is not yet known whether the storms with high continuum (with or without 
storm bursts) that begin after flares differ intrinsically from the ordinary (‘Type I) 
storms. 

Flare data have been taken from the Quarterly Bulletin, and reports issued by 
C.R.P.L. (National Bureau of Standards, Boulder, Colorado), the Fraunhoffer 
Institute (Freiburg-im-Bresgau) and the National Standards Laboratory, 
C.S.1.R.0. (Sydney). 

When it is definitely known that flare observations cited were not in progress 
during an event, a note to this effect is inserted in the Flare column; however, 
the complete absence of an entry in the Flare column could mean either that no 
observations were in progress or that no flare was observed. 

The S.1.D.’s listed have been taken from the C.R.P.L. reports. In the same 
column are occasionally listed the starting time, duration and strength of polar 
blackouts, the letters P.B. being used to distinguish them from S.I.D.’s. The 
data on polar blackouts are those of Reid and Leinbach (8), with additional 
recent information from Arnoldy, Hoffman and Winckler (9). 

The heliographic distribution of flares associated with the events of Table I 
is shown in Fig. 2. It is seen that Type II-associated flares are distributed in 
longitude across the whole disk; the gradual falling off towards each limb could 
be due wholly or partly to the fact that flares are more difficult to see near the limb 
than near the centre. The storm-associated flares are also widely distributed, 
but 20 out of the 23 fall between longitudes 0 and + 60°. 

The distribution in latitude of flares associated with Type II bursts and with 
storms both show a marked excess (about 2:1) of events North of the Equator. 
An examination of data published by Thompson and Maxwell (10) reveals a 
similar trend (18 northerly to 9 southerly) in the case of storm-associated flares 
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Fic. 2.—The heliographic distribution of flares associated with the major radio events, and 
histograms showing their east-zest and north-south distribution. 


recorded at Fort Davies, Texas, during different daily hours of observation 
within a 27-month period centred on the year 1958. Corresponding data for the 
Type II-associated flares (11) (recorded between 1956 October and 1957 
September) also show a slight, though scarcely significant, northerly increase. 

Inspection of Table I clearly shows the general increase of activity following 
the minimum of solar activity in 1954. It is interesting to note that during the 
first part of the sunspot cycle (1954-1958-5) storms with continuum were 
detected on the spectral records on only four occasions, while 75 Type II bursts 
were recorded; between 1958-5 and 1960, on the other hand, there appeared 
25 storms to 66 Type II bursts. This result indicates that the continuum storm 
may be a feature of the late phase of the sunspot cycle, perhaps requiring a 
well-developed corona for its generation. 


The observations and analysis of Table I have involved the work of different 
members of this Laboratory at different times; of the numerous people who have 
contributed, the authors wish particularly to acknowledge the major contributions 
of Mr J. D. Murray, Dr J. A. Roberts and Mr J. Joisce. 


CSLRG., 
Radiophysics Laboratory, 
Sydney, Australia: 


1961 March 3. 
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OF COSMIC RAY INTENSITY DURING THE PERIOD 
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Summary 


An unusual sequence of cosmic ray variations, observed in London in 
connection with the high level of solar and geophysical activity prevailing 
during the period 1960 November 10-16, is described. ‘The results are 
compared with cosmic ray observations made elsewhere and with other 
relevant geophysical phenomena. Some of the features encountered are 
unique in the type of event considered, and the implications of these results 
are briefly discussed. 





1. Introduction—The cosmic ray intensity observed deep in the Earth’s 
atmosphere was subject to a series of unusual changes during the period 1960 
November 10-16. Additional cosmic radiation accelerated at the Sun to energies 
high enough to be detectable at sea level was observed twice during this period, 
even though increases of this kind are relatively rare. Moreover, at least two 
Forbush decreases exhibiting a high degree of anisotropy occurred on November 12 
and 13. In addition to these cosmic ray observations, the associated series of solar 
and geophysical phenomena were observed by a wide range of experimental 
techniques. This circumstance, taken together with the remarkable super- 
position of a variety of effects which took place during this period, provides a 
valuable opportunity of deriving a large amount of information about the solar, 
interplanetary, and geophysical processes which were operating at that time. 
It is intended in this note to report the cosmic ray intensity variations observed 
with the different types of detectors located at Imperial College, and to make a 
preliminary attempt at integrating these results with similar observations made 
elsewhere as well as with the associated solar and geophysical phenomena. 

2. Solar and geophysical data.—During the week under consideration an 
exceptionally active region (McMath plage region 5925), associated with a 
complex sunspot group was passing across the visible disk of the Sun at a 
latitude of 26° to 29°N, with central meridian passage on November 12*. 
Flares of importance 3 or 3+ occurred in this region on the roth, 12th, and 15th, 
while one of class 2+ was seenonthe 14th. The flares of class 3 + of November 12 
and 15 were responsible for the additional cosmic radiation detected by sea-level 
neutron monitors at high latitude on these dates (see Section 3). A flare which 


* We are indebted to Professor M. A. Ellison who provided us with detailed records of the flares 
mentioned. 
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occurred on the 11th was missed by the visual observation network but was 
detected by its radio emission. 

In the case of the flare on November 10, maximum light output was reached 
at 10" 21™ and the associated S.I.D. showed that the ultra-violet or X-radiation 
produced was of moderate intensity. The radio emission from the flare, recorded 
over a wide band of frequencies, was of moderate intensity (~ 3 x 1071" w.m? (c/s) 
min at 200 Mc/s and 12 x 107! w.m? (c/s) min at 3000 Mc/s). 

At about 03"15™ on November 11 a flare of major importance probably 
occurred but no reports of visual observations of this event are available. 
However, strong radio emission of over 7 x 10~!’ w.m*(c/s)min at 200 Mc/s 
was reported as reaching peak intensity at 03!15™. ‘The radio emission was 
particularly strong at frequencies below 1000 Mc/s and, since strong emission 
below 500 Mc/s appears to be a distinctive feature of flares which subsequently 
produce geomagnetic storms (Sinno and Hakura 1958), it appears likely that a 
plasma cloud ejected at the time of this flare was responsible for the severe 
magnetic storm with a sudden commencement (S.C.) at 13" 48™ on November 12. 
The severity of the S.I.D. which occurred between 03" 17™ and 05" 00™ showed 
that the flare was also a strong source of ultra-violet or X-radiation. 

A flare of major importance (3 + ) occurred again in this active region (27° N, 
2 W) at 13"21™ on November 12 and the visible light from it rapidly reached 
a maximum between 13%28™ and 13%31™ (flash point). ‘The emission then 
decayed to a rather small value till 14" 20™ when a second weak outburst occurred. 
By 16" 00" the visible light had reached a low intensity even though the flare 
is listed as observable until 19"22". ‘The radio emission from the flare followed 
a very similar pattern with the observations at 3000 Mc/s showing a broad 
maximum around 13'28™ to 13"50™ and subsequent decay with a small 
secondary peak at 14"20™. The intensity maximum of 107! w.m?(c/s) 
reached at this frequency is one of the highest recorded during the current solar 
cycle. ‘This feature appears to be significant in view of the acceleration of 
particles to cosmic ray energies (several GeV*) which accompanied the flare. 
At frequencies below 1000 Mc/s the emission was not particularly strong, but 
it is worthy of note that the radiation was strongly (35 per cent) polarized during 
the later stages of the event. Radio bursts of Types II and IV were recorded at 
Sydney in association with this flare. ‘The ionospheric disturbance and the 
sudden enhancement of atmospherics (S.E.A.) which accompanied the solar 
flare were especially strong. It is probable that the plasma cloud from this flare 
was responsible for the magnetic storm with a S.C. at 10" 23™ on November 13. 

Activity recommenced in the plage region 5925 early on November 14 when 
there occurred a 2+ flare which reached its peak light intensity at 03" 04™. 
Moderate radio emission and strong S.I.D. were associated with this flare. 
The solar stream which produced the magnetic storm with a S.C. at 13"04™ 
on the 15th probably originated in this flare. 

The final great burst of activity from this plage region occurred on November 
15, when a flare of class 3 + was observed at the heliographic coordinates 26° N, 
35 W. ‘The emission of optical radiation began at 02" 07™, reached a maximum 
at 02" 21™, and subsequently decayed to very faint intensity by 04"27™. ‘The 
second increase of intensity recorded by neutron monitors during the week under 
review was associated with this flare. The burst of radio emission which occurred 


* 1 GeV = 10% electron volts. 
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during the event rose to a maximum at 02! 17™ and was strong, especially at the 
higher frequencies. Strong S.I.D. were also recorded at about this time. 

The first magnetic storm to occur during the period under consideration was 
a particularly strong one which began with a S.C. at 13"48™ on November 12. 
There can be little doubt that this storm, the S.C. of which, by a coincidence, 
occurred only 20 minutes after the optical and radio maxima of the great flare 
which generated cosmic radiation, has added much to the complexity of the 
observed phenomena. As explained earlier, the plasma cloud which arrived at 
the Earth around 13"48™ was in all probability produced by the visually 
unobserved flare at 03" 15™ on the 11th and, on this basis, the time lag between 
the flare and S.C. is 34:5 hr. If, on the other hand, the cloud is ascribed to the 
flare of the roth, the time delay would be 51-5 hr. The progress of this magnetic 
storm can be followed by reference to Fig. 1 and Fig. 4, in which the changes of 
the horizontal intensity of the Earth’s field observed at the equatorial station 
Colombo are plotted. 

About 3} hours after the S.C., at 17"00™, the intensity of the horizontal 
magnetic field at equatorial stations had decreased below the values prevailing 
before the S.C., showing that the main phase of the storm had set in. The 
decrease proceeded continuously until about 20" 00™, by which time the decrease 
amounted to over 120y. After 20"%00™ there was a small recovery till 04" 00™ 
on November 13, subsequent to which time the decrease proceeded uninterrupted 
until the minimum of over 300y below normal was attained around og"oo™. 
The magnitude of the decrease classifies the storm as an extremely severe one. 
Recovery of the equatorial horizontal field from this storm appears to set in 
after og" oo™. 

During this period of recovery, the S.C. of a second magnetic storm occurred 
at 10"23™ on the 13th, indicating that the plasma cloud from the strong flare 
of the 12th had arrived in the vicinity of the Earth taking a time of about 21 hours 
to traverse the Sun-to-Earth distance. The effect of this storm on the horizontal 
intensity at the equator was rather weak and produced only a small trough in the 
recovery curve of the earlier major storm. This magnetic storm, however, was 
associated with a pronounced Forbush decrease in the cosmic ray intensity 
accompanied by a marked anisotropy, as will be seen in Section 4 of this paper. 
The value of H at the equator continued to rise throughout the rest of the day 
and levelled off during the 14th to a roughly constant value 50 to 100 y below 
the level obtaining before the S.C. on the 12th. The intensity continued to be 
depressed during the 14th and the first half of the 15th until 13"04™ when the 
S.C. of a magnetic storm due to the class 2+ flare of the previous day occurred. 
The delay between flare and S.C. in this case is 34 hours. 

3. Observations of the solar injections of cosmic radiation.—We shall first 
consider the increases observed by the standard neutron monitor at London 
but, as already mentioned, the occurrence of a magnetic storm a short while after 
the flare on the 12th which accelerated the cosmic rays adds to the complexity 
of the observed phenomena and it is evident that the data from any one station 
can provide only a part of the complete picture of the operative factors, and 
caution should be exercised in the interpretation of the observed features. For 
this reason we shall compare our results with those from a number of stations in 
Europe which provide a useful latitude coverage over a limited range of longitudes 
before proceeding to discuss the possible origins of the variations (see Fig. 1). 
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About 30 minutes after optical and radio maxima of the major flare of the 12th, 
a significant increase of intensity was detected by the neutron monitor at London 
indicating the arrival of solar particles at the Earth. From 14500™ onwards 
the intensity showed a slow growth, which took a rather long time in reaching 
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a broad maximum of about 13 per cent above the pre-flare level around 16" 30™. 
The intensity then started a decline which continued until 17"15™ when a 
further increase of a roughly linear character commenced. It will be seen from 
Fig. 1 that the rapid diminution in H at Colombo also started at about this time. 
The slow increase continued until 19" 00™, by which time a value of 28 per cent 
above normal was being recorded. At 19"00™+ 1™ a sudden and sharp increase 
of intensity occurred with a maximum intensity of 1-6 times normal at 19"03™ + 1™. 
This sudden rise occurred approximately 30 minutes before the beginning of a 
Forbush decrease recorded by the meson telescope at sea level in London. 
Between 19"00™ and 20"15™, the neutron intensity showed a double hump 
structure as can be seen in Fig. 1. The final decline of the solar cosmic rays 
observed at London began at 20" 00™ and continued in a roughly exponential 
decay on which were superposed a number of peaks that were about 75 minutes 
apart. This series of pulses in the flux of solar cosmic rays is of a highly unusual 
character and is of exceptional interest. It may be noted that the peaks 
correspond roughly in time to a series of pulsating decreases which occurred in 
the equatorial H. By o5"00™ on November 13 the intensity of the solar cosmic 
rays had fallen below the minimum observable level at London. At no time 
during the observation of solar particles by the neutron monitor was a significant 
increase detected by the meson telescope at sea level. The production of solar 
cosmic rays, then, was confined to the low energy end of the spectrum (below 
5 GeV) during this event. 

We have included the increases observed by the neutron monitors at 
Uppsala, Leeds, Lindau, Praha, and Jungfraujoch in Fig. 1 for comparison 
with the London data. The results from Amsterdam are almost identical with 
those from London and are not shown here. The intensities are expressed as 
percentages of the mean intensity during the first half of the 12th. Whilst the 
broad pattern of variations described earlier, consisting of an increase due to 
solar cosmic rays superposed on a Forbush decrease, is evident at all these 
stations, the following dissimilarities cannot be ignored. During the initial 
phase of the increase between 14" 00™ and 16"00™ the maximum was reached 
earlier at lower latitudes than at higher ones. Also, Uppsala was observing an 
increase about twice as large as the one at Leeds. Between 17" 00™ and 19" 00™ 
the slow growth of the intensity of solar cosmic rays was most pronounced at 
London, Amsterdam, and Lindau. ‘The intensity at Uppsala was actually 
decreasing during this period, and because of this Leeds, which has a higher 
geomagnetic threshold rigidity than Uppsala at normal times, was observing 
nearly the same intensity as Uppsala just before 19"00™. ‘The double hump 
between 19"00™ and 20"o0o™ in the intensity curve occurred at most stations, 
but the initial increase at 19" 00™ was not so rapid at Uppsala and Leeds as at 
London. Finally, the pulses on the decay portion of the curves after 20 00™ 
were also most marked at London, Amsterdam, and Lindau, whilst they were 
not perceptible at Uppsala and Leeds. Moreover, the increases measured by 
the monitors at Uppsala and Leeds were still approximately equal during this 
decay phase. 

It is known that the response of a neutron monitor at sea level to the primary 
radiation falls off rather rapidly at primary rigidities less than 1 GeV*, and in 


* The rigidity is here defined as pc/Ze where p is the momentum and Z the atomic number of 
the particle. pc is expressed in GeV so that the rigidity is in units of 10° volts. 
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order to study the behaviour of the flux of solar particles at the low energy end 
it is necessary to examine the cosmic noise absorption (C.N.A.) measured at 
high geomagnetic latitudes. ‘The results from the stations Longyearbean 
(A=74°:2)*, Kiruna (A=65°-7) and Trondheimt (A=62°:5) relating to the 
event of November 12 are shown in Fig. 2 where the equivalent daylight C.N.A. 
is plotted against time, the values obtained during the night being increased by 
a factor of 3-3 in order to convert them to the equivalent daylight values. It is 
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worth pointing out that the C.N.A. is sensitive to the incidence of cosmic rays 
of low energy as well as to auroral particles of similar energies in the upper 
atmosphere, and this should be taken into account when attempting to interpret 
the results. 

It will be seen from the diagram that the solar particles of low energy first 
arrived at the Earth at about 14" 30™, that is to say, some 40 minutes later than 
the arrival time of the higher energy flux detected by neutron monitors at near-by 
stations. ‘The growth of intensity of these low energy particles was also much 
slower than that of the higher energy flux, the maximum intensity being reached 


* ) is the effective geomagnetic latitude (see Quenby and Webber 1959). 
+ We are indebted to Dr K. W. Eriksen for the data from Longyearbean and Trondheim and to 
Dr B. Hultquist for those from Kiruna Geophysical Observatory. 
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only after 04"00™ on the 13th, whilst the corresponding maximum of neutron 
intensity occurred around 16-00 on the 12th. Indeed, the entire character of 
the C.N.A. variations is radically different from that seen by the neutron monitors 
and it is evident that we are now dealing with flare particles in a different and 
much lower energy range, as is shown by the much longer time taken by the 
particles to diffuse through the inner solar system. It is also clear from Fig. 2 
that the C.N.A. variations at the three stations, spread over a wide range of 
latitudes, were very similar. Attention should also be drawn to the sharp increases 
of C.N.A. which occurred over the interval 19"00™ to 21"00™ and which were 
apparently related to the sudden increase observed by neutron monitors around 
1g"00™, The magnitude of these increases in the primary flux as detected by 
C.N.A. and as observed by neutron monitors both appear to be roughly equal 
to a factor of 3. 

A consideration of the results outlined above as well as those from other 
stations distributed on a world-wide basis, including polar stations, permits us 
to establish the following broad features in relation to this event of November 12. 

The first solar particles to arrive at the Earth appeared between 13"50™ 
and 14"00™, Although impact zones in the normal sense were not identifiable 
during the initial phase of the increase, very marked longitudinal anisotropies 
in the incident flux were present during the period 14200™ to 1800™. The 
rise in intensity was steepest and the maximum increase was largest and earliest 
at stations which, when geomagnetic deflection is allowed for, were scanning 
those portions of the sky that lie between 0° and 120° due west of the Sun—Earth 
line. In the case of stations viewing the remaining part of the sky, the growth 
was slower and the maximum intensities of solar particles recorded were only 
about half as large as those in the first group of stations. The stations in Europe 
were at this time recording solar particles arriving from directions away from the 
Sun and were therefore recording the lower intensities from these directions. 
The ratio of the intensities measured at Uppsala and London as well as the 
ratio of Leeds to London were increasing during the first two hours of the flare 
increase, and this shows that the lower energy part of the solar particles was 
building up more slowly than the higher energy portions in these directions. 

By 16"00™ to 17"00™, however, the flux of solar particles appears to have 
become sufficiently isotropic to enable a spectrum determination to be made, 
and the differential rigidity spectrum of the extra radiation from the Sun at this 
time is of the form dN/dp = (po/p)” (p in units of rigidity) where pp, is a constant 
and n = 7-2 on the average, in the rigidity range 1 to 4-5 GeV. The curve of solar 
particle intensity plotted against threshold rigidity for this event is, however, 
markedly different from that for 1956 February 23, in that the curve for the 
November 12 event flattens off more sharply at lower rigidities even though the 
spectrum appears to be steeper at the higher rigidities. The equivalent value 
of n in the recent event varies from ~5 to 1 GeV to ~8 at 4:5 GeV. ‘This could 
reflect a flattening of the flare spectrum itself. We argue that this flattening 
cannot be ascribed entirely to the fact that the low energy particles had not 
arrived in comparable abundance at that time, but was due rather to a systematic 
lowering of the geomagnetic threshold rigidities appropriate to each station. 
Such a lowering of thresholds appears to be necessary to reconcile the sharp 
differences between the intensities of low energy radiation implied by neutron 
monitor and C.N.A. measurements. A consistent explanation for the flattened 


8 
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latitude curve and the differences in intensities of the low energy particles as 
implied by the two methods of detection can be obtained by assuming that the 
geomagnetic thresholds had been lowered by a filamentary ring current of 
5x 10° amp having a radius 8 to 10 times the Earth’s radius. Such a ring 
current would also lower the equatorial horizontal field by 50y. It should be 
noted that this ring current is postulated as being relative to the conditions which 
prevailed on 1956 February 23 and is to be envisaged as a quiet time current 
which existed before the onset of the main phase of the magnetic storm. 

During the period 17"00™ to 19"00™ the intensity of the solar flux as 
recorded at high latitudes appears to decrease continuously. ‘The intensity of 
the secondary radiation measured in the geomagnetic latitude range 45° to 55°, 
however, was increasing steadily during the period. This behaviour is clearly 
seen in Fig. 1 and is borne out by the data from other parts of the world. ‘This is 
precisely the character of the changes to be expected if there was, during this 
interval, a systematic and progressive further reduction in the geomagnetic 
thresholds at all stations. It can be seen from Fig. 1 that a remarkable correlation 
exists between this reduction and the decrease of H at equatorial stations during 
the main phase of the magnetic storm. ‘The lowering of threshold rigidities 
appears as a progressive shift to higher rigidities of the knee of the curve 
depicting solar particle intensity against the normal threshold rigidity (e.g. the 
curve is now almost completely flat below 1-8 GeV). It is obvious that the 
lowered thresholds produce an increase in the primary flux entering the atmosphere 
at all stations, but this increase is reflected in an increase in neutron monitor 
intensity at middle latitudes only because of the very weak response of the high 
latitude monitors to the additional flux. 

Around 19"00™ there occurred a sharp increase of the neutron intensity 
recorded at all stations, including those at high latitudes. The time of commence- 
ment varied from station to station by as much as 10 minutes, and the times 
taken to reach maximum intensity were also widely different. It is evident from 
the foregoing discussion that this increase cannot be attributed to a sudden 
further lowering of the geomagnetic thresholds which would introduce additional 
low energy radiation from a constant intensity of solar particles present outside 
the Earth’s magnetic field. Instead, there is little doubt that a sudden increase 
occurred in the flux of solar particles arriving in the vicinity of the Earth at this 
time (outside the geomagnetic field) and that it was related to a changing 
geometry of the trapping region in the interplanetary magnetic field. It could, 
of course, be connected also with the modulating mechanism which produced 
the Forbush decrease commencing at 19"30™. The lack of isotropy prevailing 
at this time prevents us from making an accurate determination of the energy 
spectrum of this additional radiation, but, broadly speaking, the absence of 
any correspondingly sudden changes in the Uppsala to London and Leeds to 
London ratios indicates that the increase was not accompanied by any marked 
change in the energy spectrum. The absence of change in these ratios also 
lends weight to the view that this increase was not of geomagnetic origin. 

Between 19"00™ and 20"o00™ the recorded intensities at many stations 
showed marked fluctuations, with a particularly sharp trough occurring at 
19" 30™, thus producing the double hump appearance already mentioned. The 
dissimilarities between stations are rather pronounced, the first hump being a 
particularly variable feature. We believe that this interval represents a transition 
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from a region of low intensities to another where the intensities were three times 
as great, this transition being due to the outward motion of the interplanetary 
trapping field past the Earth. The transition to high intensity was irregular and 
the particle flux became very anisotropic during this time. It is worth noting 
that a strong positive pulse of ~100y occurred at 19"30™ in the horizontal 
(magnetic) intensity at equatorial stations and the possibility is not excluded, 
therefore, that the trough in neutron intensity occurring at this time, as well as 
certain other anisotropic features during this time, were due to some geomagnetic 
effect causing a sharp increase in the threshold rigidities. 

By 20"00™, however, at most stations maximum intensity was reached 
together with the establishment of a high degree of isotropy in the flux of solar 
protons arriving at the Earth. At about this time the horizontal field intensity 
at equatorial stations was also passing through a shallow minimum of about 120 y 
below normal, and it seems likely that the maximum lowering of geomagnetic 
thresholds also occurred about this time. ‘The Uppsala to London ratio had by 
now been reduced to 1-6: 1 (as against 4:6: 1 at 16" 00™) and Leeds was observing 
nearly the same increase as Uppsala. In fact, by considering the data from a 
world-wide network of stations, it appears that all stations which had thresholds 
below 2-0 GeV at normal times were registering very nearly the same intensity. 
It is estimated that a simple filamentary ring current of radius 6 to 8 Earth radii 
and of strength 10’ amp would produce the necessary lowering of geomagnetic 
thresholds which could explain the above results. Such a current would also 
reduce the equatorial H by 1ooy below that prevailing between 16"00™ and 
17800™,. ‘This estimate is evidently in fair agreement with the observations 
of H at Colombo. 

During the decay phase which set in after 20%00™ the solar particles were 
arriving uniformly from all directions, as can be seen by the close agreement 
between the increases observed by the world-wide network of stations. 
Derivation of a mathematical form for the decay of the radiation is complicated 
by further changes in magnetic thresholds, but it can nevertheless be stated that 
the highest energy particles disappeared most quickly thus producing an 
increasingly steep spectrum of the solar flux at the Earth. The exponent n of 
the differential spectrum was 8-5 +0-3 in the range 1 to 4 GeV by 02"00™ on 
the 13th. 

Additional evidence for lowered geomagnetic thresholds at various stations 
during this event is provided by the C.N.A. measurements. Considering, for 
example, the intensities at 18"00™, we note that if the threshold at Uppsala at 
this time was the normal value of 1-2 GeV as determined by the Earth’s internal 
field, the additional flux of flare particles at the top of the atmosphere could be 
deduced, using specific yield functions, to be about 100 times the normal 
galactic flux above this energy. If we now extrapolate the spectrum (with exponent 

= 7:2) derived from cosmicray observations in the range 1 to 4 GeV down to lower 
energies, we can predict that the C.N.A. at Trondheim with a normal threshold 
of o-7 GeV should be 4dB corresponding to an intensity of flare particles 
1000 times the normal galactic flux. As a matter of fact the observed absorption 
is 12 dB, indicating that the estimated flux is too low by a factor of 10 or more. 
Indeed, the estimated value of 4 dB is probably already too high because the 
lower energy radiation, as measured by the C.N.A., was still increasing at this 
time. 





8* 
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If, however, we assume that the threshold at Uppsala was 0-7 GeV the 
corresponding flux of solar particles at the top of the atmosphere is 1000 times 
the galactic flux and we would now expect to observe an absorption of 12 dB 
at stations having a threshold of 0-5 GeV or less. It appears, then, that the 
cosmic ray and cosmic noise absorption measurements can be reconciled only 
if we assume a systematic lowering of the geomagnetic thresholds appropriate 
to each station. ‘Taking this in conjunction with the cosmic ray evidence it is 
apparent that both before and particularly after the main phase each observatory 
was admitting particles having energies very much lower than could conceivably 
have been allowed on the basis of thresholds determined solely by the internal 
magnetic field of the Earth. 

Mention should also be made of the pulses of neutron intensity which were 
superposed on the exponential decay of the radiation after 20" 00™ on the 12th. 
These pulses were the strongest around the latitude of London but were not 
observed at geomagnetic latitudes above 55°. ‘They occurred at the same 
universal time at all stations which were in a suitable position to observe them 
and were centred at 21"10™, 22!40™ and oo'30™U.T. As the latitude 
dependence of the magnitude of these sharp increases is closely similar to that 
encountered earlier for the increases during the period 17"00™ to 1g9"00™, 
when a lowering of geomagnetic thresholds occurred, and because of the 
correspondence in time between these pulses of solar cosmic rays and the 
pulsating decreases observed in equatorial H, it seems very likely that these 
enhancements of intensity were caused by sudden and short-lived decreases in 
the geomagnetic thresholds which recurred at intervals of about 1} hours. 

In contrast to the great complexity of the flare outburst of November 12, 
the second flare-associated increase of neutron intensity which followed the 
3+ flare of o2"21™ on the 15th showed all the characteristics that have come 
to be regarded as typical of such events. 

At London, the neutron intensity rose gradually from o2"40™ onwards, 
reaching a maximum of about 15 per cent above the pre-flare level by about 
03" 30™ and thereafter decaying away in a nearly exponential manner. ‘These 
results, together with those from other stations in Europe, are shown in Fig. 2 
were it can be seen that similar features were observed at other stations as well. 
However, these stations were not lying within impact zones for this event. The 
rapid and large increases characteristic of impact zones were seen by stations 
receiving primary particles arriving from asymptotic directions lying between 40° 
and 100° to the west of the Earth-Sun line. ‘The oscillations reported by 
Steljes, Carmichael and McCracken (1961) were observed during the transition 
from this initial phase to the isotropic phase. The oscillations, however, appear 
to have been even more highly localized than the impact zones themselves. 

Almost complete isotropy of the incoming particle radiation had set in between 
04" 30™ and o5"00™ and about this time the exponential decay of the flare flux 
also commenced. A characteristic time ft, for the decay could be deduced, which 
varied systematically from station to station (for example, ty)= 4-0 hr at Uppsala, 
ty =3°5 hr at London, ¢,=1-8 hr at Jungfraujoch). ‘The characteristic time was, 
in fact, rigidity dependent with a form t)(p)=4:5p~”, where ¢, is in hours and 
pin GeV. ‘The rigidity spectrum at o5"00™ had the same exponent as that at 
17"00™ to 1800™ during the earlier increase. Even the flattening-off of the 
solar particle intensity versus normal threshold curve below about 1-5 GeV, 
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observed early in the November 12 event, was reproduced. ‘This last fact is 
easily understandable because the Earth’s magnetic field had not recovered 
completely from the storms of the 12th and 13th and it is likely that the 
geomagnetic thresholds were still below normal. 

The increase of C.N.A. remained below 2 dB until 06"00™, at which time 
the intensities recorded by neutron monitors had already been decaying from 
their maxima for about an hour. The intensities of low energy particles implied 
by the C.N.A. measurements indicate that the p~* spectrum of the solar 
particles, as measured at rigidities above 1 GeV, must have flattened out below 
0:8 GeV. The flux of low energy particles in fact continued to grow throughout 
the major part of the 15th and produced a maximum daylight equivalent C.N.A. 
of about 50 dB by 20" 00™, at which time the neutron monitors at high latitudes 
were showing increases of only a few per cent above the galactic flux. The solar 
particle flux as implied by neutron monitor measurements and the C.N.A. 
measurements are consistent provided we assume lowered thresholds are still 
operative at this time. We estimate that the flux at this time had reached an 
intensity of ~7x 104 per cm? sec~! above an energy of about 10 MeV. The 
value of the maximum C.N.A. suggests that the flare of the 15th was much more 
copious in its production of low energy particles than the earlier one, when 
the maximum absorption was about 30dB corresponding to an intensity 
~ 2°5 x 104/cm?sec—1. 

4. The Forbush decreases.—In addition to providing these spectacular examples 
of the acceleration of particles to cosmic ray energies during solar flares, the high 
level of solar activity prevailing during this eventful week was also responsible 
for the production of a series of Forbush decreases which followed one another 
in rapid succession. At high latitudes, these decreases were best observed by 
meson detectors which were too insensitive to detect the large influx of solar 
particles which overlapped the Forbush decreases on both occasions. ‘The 
detectors at London consisted of standard cubical telescopes operated at sea level 
and two semi-cubical scintillator telescopes having a high counting rate 
(50000/hr) operated underground at a depth of 60 metres water equivalent 
(m.w.e.)*. The data from these recorders and the neutron monitor are plotted 
in Fig. 4, which also shows the changes of the horizontal magnetic field at Colombo. 

It can be seen from Figs. 1 and 4 that the first decrease of the series commenced 
at 19230™+15™, nearly 6 hours after the S.C. and about 2} hours after the 
beginning of the main phase of the magnetic storm accompanying the decrease. 
It would appear that the commencement and the initial stages of this Forbush 
event, within the limits of experimental accuracy, were simultaneous on a 
world-wide scale. A maximum depression of about 4 per cent was reached within 
an hour or two of the commencement by the counter telescope at London, and 
thereafter a slow recovery to normal started. ‘The underground monitor also 
followed a pattern of variations similar to that of the telescope at sea level, but 
the maximum decrease was 1:'2+0°3 per cent. The rapidity with which the 
decrease took place, particularly in the high energy radiation measured by the 
underground telescope, is remarkable, and this determines the strength of the 
screening magnetic field that is required to move past the Earth as being of the 
order of 100y if we assume a mean rigidity of 200 GeV for the primary particles 


* The absorbing material over the recorders has the same effect as a water absorber of 60m 
thickness. 
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which are monitored at a depth of 60 m.w.e. Comparison of our results with 
those from other stations shows that the recovery from the first decrease was not 
associated with any marked anisotropies. 

This first decrease cannot be readily identified in the intensities registered by 
neutron monitors at latitudes over 50° because of the prevailing flare-associated 
increases, but it is possible to deduce that the decreases in these monitors must 
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have been about 6 to 8 per cent. The flux of solar particles had decayed to less 
than a few per cent even at high latitudes by 12" 00™ on the 13th, and between 
this time and 02! oo™ on the 15th the intensity variations recorded by the neutron 
monitor and the cubical telescope at London were roughly parallel. 

A short while after 10" 30™ on the 13th, a second Forbush decrease occurred 
which depressed the intensities further in all three recorders at London. Unlike 
the first decrease, the second one began at about the same time as the S.C. of the 
weak magnetic storm which was probably caused by the arrival of the plasma 
cloud from the flare of the 12th. The time taken to reach the minimum between 
21%00™ and 23"00™ was much longer for this event, and during the course of 
the decrease exceptionally strong anisotropies were developing. The solar 
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particles that were still arriving at high latitudes from the event on the 12th were 
completely removed at the onset of this event although no corresponding depletion 
of the low energy particles as measured by the C.N.A. was noted at this time. 

The total decrease resulting from the two Forbush events amounted in London 
to 13 per cent in the neutron intensity, 8-5 per cent in the meson intensity at sea 
level, and 1-4 per cent in the meson intensity at a depth of 60 m.w.e. The ratio 
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of 0-15 +0°03 for the decrease measured underground to that in meson intensity 
at sea level would imply, on the basis of specific yield functions available at 
present, that the fractional modulation of the primary intensity during this event 
could not have been steeper than a power of —1 in rigidity (or energy) over the 
range 20 to 200 GeV. In fact the average modulation over this range fits the 
results closely when represented by a function of the form 

dP_ Py (dy 

aa ~ Pa) 
where dJ,/dP is the unmodulated spectrum and Py is a constant. 

During the early hours of the 14th the recovery of the meson intensity had 
commenced and the period of recovery was attended by continued strong 
anisotropies. Large anisotropies were also observed by the underground 
detector. There is some evidence that the axis of the anisotropy was such 
that the greatest decrease occurred at a local time of 03" 00™, and this direction 
of anisotropy was maintained until the anisotropy disappeared on the 15th. 

These features are identical for the very high energy cosmic rays observed by 
the underground detector as well as the lower energy particles observed by 
neutron monitors at sea level. 
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It is of interest that the direction from which solar particles arrived in direct 
paths from the Sun in this second flare outburst was roughly perpendicular to 
the axis of anisotropy of the Forbush decrease which was still in progress at this 
time. 

A more detailed analysis of various aspects of these Forbush decreases, 
including a comparison with models suggested as explanations of these events, 
is at present under way. 

5. Conclusions.—The remarkable sequence of events described above gives 
promise of much valuable information on the conditions existing in interplanetary 
space and in the outer regions of the Earth’s magnetic field at this time. 
Meanwhile, it is convenient to summarize some aspects of the flare increases 
which have become clearly evident for the first time in the course of this event 
and which can provide quantitative information on the factors controlling the 
propagation of solar cosmic rays between the Sun and the Earth. 

(i) The important effects on the Earth’s magnetic field (Webber 1961) and 
possible modifications in its structure were amply demonstrated by such results 
as the apparent flattening of the flare energy spectra at low energies and the 
apparent intensity variations which did not necessarily reflect the conditions 
existing outside the influence of the Earth’s field. 

(ii) ‘The first event in particular showed that trapping regions with well-defined 
boundaries may exist in the inner solar system which may be filled to much higher 
intensities with solar particles than the surroundings. There is quite probably a 
relationship between these trapping regions and the mechanism producing 
Forbush decreases. ‘There is also the possibility that the trapping regions may 
be linked directly to the flare regions. 

(iii) ‘The intensity-time characteristics of both the initial and the decay 
stages of these events, as well as the degree of anisotropy, which was ~ 1 per cent 
during the decay phases of both events, indicate that scattering and diffusion are 
important aspects of the interplanetary magnetic field between the Sun and the 
Earth as well as outside the Earth’s orbit. 

(iv) The preferential arrival of solar particles from directions lying well to 
the west of the Sun was encountered at the beginning of both events, as also during 
the event of 1960 May 4 (McCracken and Palmeira 1960), and this suggests the 
action of a deflecting influence which changes the apparent position of the 
“source”. ‘This preferential direction may be closely related to anisotropies 
existing in the galactic radiation at these times. 

(v) The flux of sofar particles at low energies may reach a level 10° times the 
normal cosmic ray flux, and the corresponding energy density about 1000 times 
normal or 1000 x 107!” erg/cm?. 
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Summary 


Part of the galactic plane has been surveyed at 408 Mc/s using the 250 ft 
radio telescope at Jodrell Bank. The results are compared with the Sydney 
85 Mc/s survey and the Leiden 1390 Mc/s survey in order to group the 
discrete sources of radio emisson according to their spectral index. Of fifty- 
one sources, thirty-six are found to have a spectral index consistent with the 
assumption that they radiate thermally, and fifteen, including the source at 
the galactic centre, have a non-thermal spectral index. Most of the sources 
lie within 1/4° of the new galactic plane, and are thought to be at distances of 
several kiloparsecs. Some of the thermal sources are a few degrees from the 
galactic plane, and probably less than one kiloparsec from the Sun. The 
continuum background emission from the galactic disk is discussed in a 
second paper. 





1. Introduction.—The radio emission from the Galaxy has been divided into 
three major components. First, the corona (Baldwin 1955), which has a very 
extended distribution and is responsible for the greater part of the emission at 
metre wavelengths. Second, a disk component (Mills 1959), which is closely 
confined to the galactic plane, and is roughly co-extensive with the visible parts 
of the Milky Way. Third, a population of discrete sources of radio emission 
(Hanbury Brown and Hazard 1953; Mills 1959) which lie close to the ridge line 
of the disk. A large number of these discrete sources are thermal, having a 
spectrum characteristic of a hot gas. Some of these sources have been identified 
with emission nebulae (Westerhout 1958). ‘The spectrum of the remainder 
suggests that the radio emission is produced by a non-thermal mechanism. The 
three known supernova remnants in the Galaxy are radio sources of this kind, 
and the suggestion has been made that all non-thermal galactic sources are super- 
nova remnants (Mills 1959; Hanbury Brown 1954). 

In this paper the results of a 408 Mc/s survey of part of the galactic plane 
are presented. The survey is used in conjunction with the Sydney 85 Mc/s 
survey (Hill, Slee and Mills 1958) and the Leiden 1390 Mc/s survey (Westerhout 
1958) to investigate the spectrum and distribution of the discrete sources in the 
area common to the three surveys. ‘The continuum background emission from 
the galactic disk will be discussed in a second paper. 


* Now at C.S.I.R.O., Sydney. 
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2. Observations.—In May 1959 the region of the sky bounded by b= +6° 
and /¥=350° to 56° was surveyed at 408 Mc/s using the 250ft radio telescope 
at Jodrell Bank (Lovell 1957)*. ‘The equipment was substantially the same as 
that used in previous surveys (e.g. Large, Mathewson and Haslam 1959). ‘The 
aerial beam had a width to half power points of 40 x 54 minutes of arc. The 
beam solid angle as defined by Seeger, Westerhout and van de Hulst (1956) was 
measured to be 0-70 square degrees. ‘The receiver was stabilized by a Dicke 
system; the overall noise figure was 6dB and bandwidth 4 Mc/s. A daily cali- 
bration was made by observing the Cygnus A (19N4A) radio source, taking the 
flux density of that source to be 4210 x 10-7 w.m.~*(c/s)~! (Seeger, Westerhout 
and Conway 1957). ‘The region was surveyed mainly by driving the telescope 
in right ascension at 1 deg/min at constant declination. Fig. 1 is typical of these 
scans, which were made at 20 minutes of arc intervals. A few scans along the 
galactic plane were made by driving the telescope in galactic longitude, at constant 
latitude. An example is shown in Fig. 2. 
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Fic. 1.—A typical scan in right ascension at constant declination, 8= —11° 50’. 


Fig. 3 (a) shows the 408 Mc/s radio isophotes drawn in units of 43 °K brightness 
temperature, measured above a selected reference region (at «= 19" 20™, 5=0°). 
The brightness temperature is calculated from the chart deflection by the method 
described in a previous paper (Mathewson, Large and Haslam 1960). 

The position of the beam has been computed from the azimuth and elevation 
dial readings of the telescope. Corrections have been added for refraction, and 
for the deviation of the radio axis from the nominal axis of the paraboloid. ‘These 
corrections were determined experimentally by observing strong radio sources 
at various elevations. The probable error of position over the whole survey is 
+5 minutes of arc. The error in the quoted position of weak or extended 
sources may be greater than this value. 


3. The list of discrete sources.—Fig. 2 shows a 408 Mc/s radio profile obtained 
by scanning the telescope along the path close to the galactic plane. ‘The discrete 
sources are clearly seen above the background. In measuring the sources a 


* The symbols 41 and /@ are used to indicate the new galactic latitude and longitude (Trans. 
I.A.U., 1958). 
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smooth curve has been drawn through the most marked minima of such profiles 
to represent the background component. A similar procedure has been applied 
to the 85 Mc/s and 1390 Mc/s surveys to obtain a consistent definition of a source, 
and in order to compare the source flux densities at the three frequencies. 

Table I is a list of all the discrete sources recognized within the limits of the 
408 Mc/ssurvey. They have been given reference numbers in order of increasing 
longitude. Subsequent columns of the table indicate, where appropriate, the 
source number assigned by Westerhout (1958), by Wilson and Bolton (1960), 
by Mills, Slee and Hill (1958 and 1960), by the third Cambridge survey (Edge 
et al. 1959) and by the I.A.U. 

Most of the sources appearing on Hill, Slee and Mills’ 85 Mc/s survey and on 
Westerhout’s 1390 Mc/s survey have been detected at 408 Mc/s, and the table 
lists the position measured from the 408 Mc/s results. Occasionally, when a 
source is barely detected at 408 Mc/s, the position has been measured from Hill, 
Slee and Mills’ survey. The positions of these sources appear in italics in the 
table. The positions of the sources at each frequency are marked on the key 
map (Fig. 3()). It will be seen that the positional agreement of the 408 Mc/s 
survey with either the 85 Mc/s or the 1390 Mc/s is good but that frequently 
sources do not appear on all three surveys. The next columns of the table 
show 408 Mc/s flux density and the measured brightness temperature (7,) of 
each source at each frequency. In the case of the 1390 Mc/s results, the 
temperature quoted has been obtained by convolving the published 1390 Mc/s 
survey (made with a 35 minutes of arc beam) with a 32 minutes of arc gaussian 
function to render the effective resolution equal to that of the other two surveys. 
The ratios of these three brightness temperatures may then be used to obtain a 
spectral index of brightness temperature for a source, and thence a spectral index 
of flux density, which is tabulated. A study of these spectral indices (in Section 4) 
has led to a division of the thermal sources into four groups with different 
spectral properties, and has also distinguished clearly the non-thermal sources. 

4. The continuum spectrum of the discrete sources.—The spectral index of flux 
density for a thermal radio source may range from +2 for optically opaque 
regions to —o-1 for optically thin regions. In the sky, thermal sources are fre- 
quently associated with emission nebulae (H 11 regions) where the ionization of 
the hydrogen is maintained by the intense ultra-violet radiation of O and B stars. 
In such nebulae, the electron temperature is thought to be usually about 104 °K 
(Aller 1956). If other mechanisms of ionization are present, such as collisional 
excitation, then the electron temperature can be much higher than this (for 
example in NGC 6992, 6995 it is about 4x10*°K (Oosterbrook 1958)). 
Since a thermal source absorbs incident radiation, the apparent brightness 
temperature (7J,) depends on the background brightness temperature 7',, and 
is given by 

T,= (T.- Ty)(1 a 

where 7, is the electron temperature of the gas, and 7 is its optical depth (Mills, 
Little and Sheridan 1956). ‘This shows that if 7, is comparable with 7, then 
the excess temperature of the thermal source is small, and indeed it appears in 
absorption (as a negative source) if 7, is greater than 7,. Non-thermal sources 
are assumed not to absorb radiation, and thus to appear in emission against any 
background temperature. ‘The spectral index of galactic non-thermal sources 
is observed to range between o and —1 (Whitfield 1959). 
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Fic. 3.—408 Mc/s isophotes (a), and key diagram (6). 


The intensity unit is 43 °K brightness temperature. 
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The tabulated source positions are marked with dots on the isophotes. 
non-thermal are indicated by closed circles. 
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sources and closed circles indicate non-thermal sources. The three sources discussed in Section 4 (iti) which may be thermal or 
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Fourteen sources in Table I have spectral indices lying between —o-2 and 
—o-7, and are therefore classed as non-thermal (the source 17S2A at the galactic 
centre is considered separately). Some of the sources are detected at 85 Mc/s 
only, or at 85 Mc/s and 408Mc/s only. ‘This seems to be because at the 
higher frequencies weak non-thermal sources are not easily distinguished from 
the thermal emission of the Galaxy. Where measurements were possible on 
the three frequencies it is seen from Table I that consistent values of spectral 
index were obtained. None of these non-thermal sources has been optically 
identified. 

The remaining thirty-six discrete sources in Table I have spectral indices 
between —o-1 and +2:0 and are consequently classed as thermal. It is very 
likely that most of these sources are H II regions, and a number of them have been 
so identified by various authors. 

Differences of spectral index observed in the present work are discussed by 
splitting the thermal sources into four groups, as follows: 


(i) T'wenty-nine of the sources have spectral indices close to zero as measured 
between 1390 Mc/s and 408 Mc/s, and appear weakly in emission or absorption 
at 85 Mc/s, according to whether the background temperature at that frequency 
is greater or less than 10‘°K. Many of them are not detected at 85 Mc/s. 
They are probably all H11 regions with typical electron temperatures of 10* °K. 


(ii) Two of the thermal sources (Nos. 24 and 43) appear in absorption against 
an 85 Mc/s background of less than 10* °K, and thus appear to be H 11 regions of 
abnormally low electron temperatures (less than 7ooo°K). A similar result 
was reported by Mills, Little and Sheridan (1956) for the source No. 2 (W22), 
but the revision of Mills’ (1958) temperature scale renders this result invalid, 
and there is no evidence that the electron temperature of W22 differs appreciably 
from 1o* °K. 


(iii) Three sources (Nos. 13, 25, 29) have spectral indices close to zero, and 
yet appear strongly in emission against a high 85 Mc/s background temperature. 
They are either non-thermal sources with very flat spectra indeed, or thermal 
optically thin regions with abnormally high electron temperatures, say 2 x 104 °K. 
It is generally thought that H11 regions do not have electron temperatures as 
high as this, but such high temperatures can be achieved by collisional excitation 
mechanisms. Observations at a much lower frequency than 85 Mc/s would be 
needed to determine the nature of these objects. They do not coincide with any 
obvious emission nebulae. 


(iv) Two sources (Nos. 21, 35) have spectral indices apparently greater than 
zero, indicating a source optically opaque at 408 Mc/s. By assuming the true 
electron temperature to be 104 °K, and calculating the optical thickness from the 
spectral index, it is possible to calculate the brightness temperature of the source. 
This is compared with the apparent brightness temperature to find the dilution 
factor in the aerial beam and hence the source diameter. In the case of source 
No. 35 (W43) this diameter was found to be 5 minutes of arc, in reasonable agree- 
ment with an interferometer measurement (Wilson and Bolton 1960). 


5. The distribution of the discrete sources in latitude.—The distribution in 
galactic latitude of the sources listed is illustrated in Fig. 4, in which the number 
of sources lying in each $° range of latitude is plotted. One histogram includes 
the sources for which the spectral index lies in the range +2 to —o-1, and the 
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other histogram the sources for which the spectral index lies between —o-1 and 
—1-o. In the following discussion, the sources in the first group are called 
thermal, and those in the second group non-thermal. It is seen from this 
diagram that the majority of the sources lie within 1/4° of the galactic plane. 
Only one non-thermal source lies more than 1° from the plane. A number of 
thermal sources are found several degrees from the plane. 
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Fic. 4.—The upper histogram shows the number of sources with spectral index greater than or 
equal to —o-t falling in each }° range of galactic latitude. These are the sources assumed to be 
thermal in the text. In the lower histograms the sources of spectral index less than —o:1 are plotted. 
These are the non-thermal sources. 


5. (a) The thermal sources.—The distribution of the thermal sources may be 
described by means of a model distribution in which approximately one third of 
the sources are evenly distributed in latitude up to + 5° from the galactic plane, 
and the remaining two thirds have a gaussian distribution with standard deviation 
0°-35, symmetrical about the galactic plane. The galactic plane which is used 
in this description of the latitude distribution of the sources is that given by b! =o, 
which is determined by the large scale distribution of neutral hydrogen in the 
Galaxy (Gum, Kerr and Westerhout 1960). ‘The observation that most of the 
thermal sources lie very close to this plane suggests that they too are a large scale 
feature of the Galaxy. The neutral hydrogen in the Galaxy lies in a disk about 
200 pc thick at the half density points. The observations of Van Tulder (1942) 
and Schmidt (1959) have shown that, in the vicinity of the Sun, O and B stars 
(which maintain the ionization in H1I regions) are slightly more concentrated to 
the galactic plane, the distribution being about 150 pc thick. High frequency 
radio observations of the thermal emission indicate that the disk of ionized 
hydrogen is about 200 pc thick, appearing in the sky as a band of radiation 1°-6 
wide (Westerhout 1958). It is therefore surprising to find that a majority of the 
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discrete thermal sources in Table I lie within 1/4° of the galactic plane. ‘Two 
factors may account for this extremely narrow distribution of the thermal sources 
in latitude. It is possible that there is an observational bias because weak sources 
are less readily detected when they lie on the gradients away from the galactic 
plane. Also, a number of the sources classed as thermal may in fact be non- 
thermal sources of extremely flat spectra. We conclude that the majority of 
these thermal sources are individual H 1 regions or fine structure in the more 
continuous thermal emission of the Galaxy, and at distances of a few kpc. 

5. (b) The non-thermal sources—The non-thermal sources fit a model 
distribution which is gaussian with standard deviation 0°-65. Since there are 
comparatively few non-thermal sources (i.e. 14) this figure is strongly influenced 
by individual sources, and if for example the source at 54 = —1°-g is omitted 
the standard deviation is nearly halved. In interpreting the distribution of the 
non-thermal sources we use arguments similar to those used for the thermal 
sources. The symmetry of the distribution about b™=o is an indication that 
the sources are within the Galaxy, and the narrowness of the distribution can 
only mean that they are at distances of several kiloparsecs. ‘They are also com- 
paratively rare objects. ‘The mean distance from the galactic plane of six radio 
sources thought to be supernova remnants (Cassiopeia A, Crab Nebula, Cygnus 
Loop, 1C443, Tycho Brahe’s and Kepler’s Supernovae) is about 120 pe, i.e. 
they are strongly concentrated to the galactic plane. ‘Thus the distribution of 
the fourteen non-thermal sources is consistent with the suggestion that they are 
distant supernovae. 

It is not worth while to proceed to a more detailed calculation until more 
information on the range and space density of type I and type II supernovae is 
available, but it is interesting to note that our non-thermal sources have a distri- 
bution very different from that of ordinary novae in the Galaxy, for the mean 
z-distance of 80 observed novae (Payne-Gaposchkin 1957) is over 500 pc. Such 
a value cannot be taken as typical of the z-distance of the galactic non-thermal 
sources because it would imply a distribution in latitude far wider than that which 
is observed. 

6. The source at the galactic centre—The present 408 Mc/s observations of 
Sagittarius A seem to be consistent with earlier higher resolution studies at 
1420 Mc/s (Davies and Jennison 1960) and 8000 Mc/s (Drake 1959), which 
show that there are several small diameter sources within }° of the galactic 
centre. One of these sources (at /1=0°-18, b4=0°-07) Drake claims is non- 
thermal. Detailed examination of the 408 Mc/s isophotes at this position shows 
no trace of a point source, and this makes it unlikely that Drake’s source has a 
non-thermal spectrum. 
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(Received 1961 June 8) 


Summary 


The 408 Mc/s Jodrell Bank survey of the galactic plane has been used in 
conjunction with the Sydney (85 Mc/s) survey and the Leiden (1390 Mc/s) 
survey to delineate the distribution of thermal and non-thermal radio emission 
from the Milky Way. ‘The results of this analysis are presented in the form of 
both isophotes and profiles at constant galactic latitudes. They extend and 
confirm a similar analysis carried out by Westerhout. Some features of the 
distributions could indicate a spiral structure of the Galaxy, but it is found 
that there are difficulties in fitting Mill’s spiral model to the parts of the Milky 
Way seen from the northern hemisphere. 





1. Introduction.—In the preceding paper (Large, Mathewson and Haslam 
1961), which will be referred to as Paper I, the results of a survey of part of 
the galactic plane at 408 Mc/s using the 250 ft radio telescope at Jodrell Bank 
were described. In that paper the isophotes were presented, together with a 
discussion of the discrete source component of the galactic radio emission. 
A study has now been made of the more extended emission from the galactic 
disk, the 408 Mc/s results being compared in detail with the Sydney 85 Mc/s 
survey (Hill, Slee and Mills 1958) and the Leiden 1390 Mc/s survey (Westerhout 
1958). The variation of spectral index with galactic latitude in regions close 
to the galactic plane has been used by Westerhout (1958) as evidence for a 
concentration of ionized hydrogen radiating thermally. Mills (1959) identifies 
certain “‘ steps” in the metre-wave radiation from the Galaxy with points 
where the line-of-sight is tangential to the spiral arms. The purpose of this 
paper is to present the results of a new spectral analysis of the galactic radiation, 
and to interpret these results in the light of the conclusions of Mills and 
Westerhout. 

2. Theory.—An inspection of surveys of the galactic plane on various 
frequencies shows that the brightness temperature varies less rapidly with 
frequency at points within 1° of the galactic plane than at points at higher 
latitudes. It is known that in other parts of the Galaxy the spectral index is 
remarkably constant (for example, see Costain 1960) as it also appears to be in 
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the Andromeda nebula (Baldwin and Costain 1961). Shain’s survey at 19-7 Mc/s 
(1957) shows a minimum in the galactic radiation in the position of the galactic 
plane, which is thought to be caused by ionized hydrogen emitting thermally 
and appearing in absorption against the very bright background temperature. 
Accordingly, we make the assumption that the variation of spectral index over 
the galactic plane is due to a variation in the proportion of the emission that is 
produced by thermal processes. While this is a reasonable assumption it is 
important to bear in mind the possibility (stressed by Mills (1959)) that some 
or all of the variation may represent a real variation in the spectrum of the 
non-thermal emission. 

Westerhout (1958) has shown that if certain other assumptions are made, 
then the apparent brightness temperature, 7, of a mixture of thermal and 
non-thermal sources is given by 


T=(T,+T,/7) (1—e~*), (1) 


where 7, is the electron temperature of the thermal material, 7 is its optical 
depth, and 7,, is the brightness temperature of the non-thermal emission 
(supposing thermal material to be absent). Westerhout applied this expression 
to his own survey and to the 85 Mc/s Sydney survey to separate the thermal 
and non-thermal contributions to the 1390 Mc/s brightness temperatures. 
A similar analysis has been carried out incorporating the 408 Mc/s results 
(Section 4). The assumptions on which this analysis is based are briefly 
discussed : 

(i) The variation of spectral index across the galactic plane is due 
principally to the emission from ionized hydrogen. 

(ii) The electron temperature of the ionized hydrogen is constant and 
equal to ro*°K. Unless the electron temperature is grossly different from 
this value, it scarcely affects the analysis. 

(iii) The sources of thermal and non-thermal emission are well mixed in 
the line-of-sight. This condition is important only when 7 ~ 1, and thus applies 
only to the 85 Mc/s results near the central regions of the Galaxy. 

(iv) The non-thermal sources do not appreciably absorb incident radiation. 

(v) The spectral index of the non-thermal radiation is unvarying with 
frequency between 1390 Mc/s and 85 Mc/s. 

In applying equation (1) to the results of the surveys at 85 Mc/s, 408 Mc/s 
and 1390 Mc/s, the optical depth 7 has been taken to be inversely proportional 
to the square of the frequency, and the non-thermal brightness temperature to 
have a spectral index «, thus: 


7 o f-? 
In principle, the results of surveys at three frequencies could be used to deduce 
distribution of three quantities, such as 7 and T,, at one frequency, and the 
value of «. The details of the analysis are described in the following sections. 
3. Calibrations—The method of observation, and details of the receiver 
will be found in Paper I. For the present purpose we are particularly concerned 
with the accuracy of the temperature scales of the three surveys used in the 


analysis, with the determination of the true zero level of brightness temperature 
and with the contribution from extragalactic radiation. 
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The spectral analysis of the discrete sources (Paper I) provides good evidence 
that the temperature scales of the three surveys are in the correct ratio; for 
the spectral index of many of these sources was found to be close to zero, an 
appropriate value for the optically thin Hu regions with which these sources 
were identified. Furthermore, the non-thermal discrete sources had spectral 
indices which were consistent as measured between 85 Mc/s and 408 Mc/s and 
408 Mc/s and 1390 Mc/s. 

Since the flux densities of the discrete sources were measured by using 
nearby regions in the galactic plane as reference points, further information is 
required to determine the absolute brightness temperatures. Westerhout (1958) 
estimates that the true brightness temperatures at 1390 Mc/s are obtained by 
adding 1-5 °K to the temperatures indicated by his published isophotes, over 
the greater part of his survey. He obtained this figure by using the goo Mc/s 
survey of Denisse et al. (1957) and assuming a suitable spectral index. The 
brightness temperatures at 408 Mc/s were measured relative to the temperature 
of a reference region (at «=19" 20™, §=0°); therefore it is necessary to 
determine the true brightness temperature of this point. In the analysis 
which follows, this temperature (Tp) of the reference region is taken to be 


Trrp=80+15 °K 


at 408 Mc/s. It was determined by first measuring the reference region 
temperature relative to cold parts of the sky at high galactic latitudes, and 
estimating the temperature at these points by using the Sydney 85 Mc/s survey. 
The result was confirmed by repeating in detail the method used by Westerhout 
in establishing the true temperatures in his survey. No zero level correction 
is applied to the 85 Mc/s results. 

Westerhout has pointed out that the 85 Mc,s brightness temperatures are 
more consistent with his and other surveys if the temperatures at 85 Mc/s are 
multiplied by a factor of 1-32. No such correction has been applied in the 
present analysis. Indeed the 408 Mc/s results seem to be more consistent 
with the published 85 Mc/s temperature scale. 

No attempt has been made to correct any of the results for the contribution 
to the brightness temperature made by integrated extragalactic radiation. 

4. Method of analysis.—On the isophotes of radio brightness temperature 
at 85 Mc/s, 408 Mc/s and 1390 Mc/s a grid of points was constructed by drawing 
the lines of constant new galactic latitude and longitude at 2° intervals. At each 
point a “‘ background temperature”’ was read from the isophote diagrams. 
The contribution from any discrete source was not included, and any local 
irregularities in the isophotes were allowed for by a ‘“‘ smoothing” of the 
isophotes in the vicinity of the point. The results were tabulated, zero level 
corrections of 1-5 °K and 80°K being added to the 1390 Mc/s and 408 Mc/s 
figures respectively. These temperatures are taken to represent the distribution 
of radiation from the large-scale features of the Galaxy. They are presented 
in Fig. 1 in the form of smooth profiles at constant galactic latitude. 

These results have been used to determine the thermal and non-thermal 
contributions to the brightness temperature at 408 Mc/s by applying the 
theory described in Section 2. In order to fit the data to equations (1) and (2) 
a table was computed (using the Manchester Mercury Computer) giving the 
brightness temperatures at the three frequencies as a function of the spectral 
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Fic. 1.—Profiles showing the longitude distribution of galactic radio emission at 85 Mc/s, 
408 Mc/s and 1390 Mc/s at five latitudes. The profiles are derived from the Sydney, Jodrell Bank 
and Leiden surveys in the manner described in the text. The vertical scale is the brightness tem- 
perature at 408 Mc/s. The temperature scale at 85 Mc/s is obtained by multiplying these numbers 
by 25, and at 1390 Mc/s by dividing by 12:5. 


index, a, of the non-thermal component. Using the table and the three 
brightness temperatures at a point, it was possible to determine the spectral 
index to the nearest tenth of a unit, and the thermal and non-thermal contribu- 
tions to the 408 Mc/s temperature to an accuracy of a few degrees. The value 
of « varied from point to point, but it was soon found that the variations in « 
were not significant in view of the experimental errors. At most points, a value 
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Fic. 2.—Profiles showing the longitude distribution of radio radiation at 408 Mc/s at five 
latitudes (labelled O) and its separation into thermal and non-thermal components (labelled T and N 
respectively). 


of « of 2-6 gave the best fit to the measured brightness temperatures. This 
value was then used at every point of the 2° grid, and the fit to equation (1) 
was within the experimental error for nearly all the points. The results of the 
analysis are shown in the form of profiles at constant latitudes of the thermal 
and non-thermal emission at 408 Mc/s (Fig. 2) and in the form of isophotes 
in Fig. 3. 

The observations at 408 Mc/s (and at 1390 Mc/s) are confined to galactic 
longitudes north of the galactic centre. In an attempt to extend the spectral 
index measurements to southern galactic longitudes, we have used the 600 Mc/s 
survey of Piddington and Trent (1956). The resolution of this survey is less 
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than that of the surveys we have considered so far, and there is reason to believe 
that the temperature scale is inconsistent with other surveys. We find that the 
600 Mc/s survey is in agreement with the other results in the northern sky 
if the brightness temperature scale is multiplied by a factor of 0-75. The 
contributions of thermal and non-thermal emission at 600 Mc/s along the 
line 6 =o were determined by fitting equation (1) to the 85 Mc/s and 600 Mc/s 
results in both the northern and southern sky. In order to carry out this 
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Fic. 3.—Isophotes showing the distribution of thermal and non-thermal radiation in the Galaxy at 
408 Mc/s. The intensity unit is 1°K of brightness temperature. 
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Fic. 4.—The upper curve is a profile of the radio emission from the Galaxy (at b!! =o) derived 
from the 600 Mc/s survey of Piddington and Trent (1956). The brightness temperatures have been 
mutliplied by 0-75. The middle and lower curves show the division of this radiation into non-thermal 
components. 


analysis, the 85 Mc/s map was convolved with a 3° gaussian function to make 
its resolution equal to that of the 600 Mc/s survey. ‘The result is shown in 
the form of a profile in Fig. 4. The detailed results shown in Figs. 2, 3 and 4 
depend on the assumptions listed in Section 2, and to some extent on a 
subjective judgment in reading the isophote diagrams and fitting the temperatures 
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to equations (1) and (2). However we are confident that the broad features of 
the results are correct. In the following discussions we assume that the thermal 
and non-thermal brightness temperatures are distributed as shown in the 
diagrams. 

5. The spiral arms of the Galaxy 

In this section we review the large-scale structure revealed in the thermal 
and non-thermal radiation from the Galaxy, and consider the possibility that 
it reflects the spiral structure. 

5a. The distribution of the thermal emission.—If the source of thermal emission 
in the Galaxy were concentrated into a spiral or circular arm structure, its 
presence would be revealed in the distribution of thermal radiation in galactic 
longitude. The thermal emission would be at a maximum in those directions 
where the line-of-sight passes through the greatest depth of material, i.e. in 
directions which are tangential to the arms of the Galaxy. This corresponds 
to a special case (isotropic emission, m =o) in the idealized model of the Galaxy 
proposed by Hanbury Brown and Hazard (1960). 

The distribution of thermal radiation shown in Fig. 3 indicates broad maxima 
in the emission at /4=25° and /41=o°, with a steady decrease in intensity for 
longitudes greater than 25°. Westerhout’s (1958) analysis led to a somewhat 
similar result as regards the maximum at /41=25°, but in the direction of the 
galactic centre he did not find a second maximum. 

We have analysed the 600 Mc/s survey of Piddington and Trent (1956) 
in conjunction with the other three surveys (see Section 4) in order to investigate 
the distribution of thermal radiation south of the galactic centre. Although the 
beam size of the 600 Mc/s survey was 3°, the peak in the thermal emission 
at /1— +25° is clearly indicated in the result (Fig. 4) and a symmetrical feature 
appears at about /!!= — 25° (i.e. 335°). The striking fall in the radio intensity 
at 600 Mc/s at this longitude is seen to be partly thermal and partly non-thermal. 
Thus what evidence there is supports the suggestion that at /™= + 25° the line- 
of-sight is tangential to an approximately circular arm of radius 3:5 kpc. 
This was suggested by Westerhout (1958), who deduced a radial distribution 
of ionized hydrogen in the Galaxy. Mills (1959) points out that the data are 
inadequate to prove the existence of this 3-5 kpc arm, and that the apparent 
structure could equally well be explained in terms of local irregularities. 

It is significant that apart from the maximum at /"= 25°, a slight rise towards 
the galactic centre and the well-known thermal emission region of Cygnus X 
(at /1=80°) there is no indication of any other large-scale structure in the 
thermal emission from the galactic disk between /7=0 and /Z=go°. In fitting 
a spiral arm model to these features, arms would be placed in radii of 3-5 kpc 
and 8 kpc in the northern sky. There is no indication of thermal emission 
from an arm of an intermediate radius. 

The isophotes of Fig. 3 show vividly the concentration of the thermal 
emission to the galactic plane, the mean width to half-intensity points being 
1°-5. This width is consistent with the theory that the source of the emission 
is a disk of ionized hydrogen concentrated to the plane of the Galaxy and 
having a thickness to half-density points of 200 pc. 

Between longitudes /¥=20° and /'=35°, the thermal emission appears to 
extend weakly towards positive galactic latitudes. This is the region where a 
continuation of the north polar spur of the Galaxy would intersect the galactic 
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plane (see Hanbury Brown, Davies and Hazard 1960). Thus there is some 
indication that the base of the spur has a thermal spectrum. A recent survey 
in Cambridge (Costain 1960) seems to indicate that at higher latitudes the 
spectrum of the spur is if anything more steeply non-thermal than the rest 
of the galactic corona. 

5b. The distribution of the non-thermal radiation.—In a recent paper, Hanbury 
Brown and Hazard (1960) have calculated the distribution of radio brightness 
which would be observed from the position of the Sun for spiral (and circular) 
arm models of our galactic system. The emission is taken to be anisotropic, the 
intensity depending on (cos @)” where n lies between o and +2 and @ is the angle 
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Fic. 5.—This diagram is reproduced from Hanbury Brown and Hazard’s recent paper (1960). 
It shows the expected distribution in longitude of galactic radio emission from a circular arm of radius 
6 kpc, the emissivity in a direction 0 to the radius being proportional to cos" @. Curves are drawn 
for n=0, 1, 1°7, 2. 


between the line-of-sight and the tangent to the spiral arm. Fig. 5, which is 
reproduced from their paper, shows the longitude distribution of radio emission 
expected from an arm of 6kpc radius for various values of n. If, as is widely 
believed, the non-thermal emission of the galaxy is synchrotron emission, then n 
would be about 1-7. This supposes perfect alignment of the magnetic field along 
the spiral arms. If random deviations occur, rather more pronounced “‘steps”’ 
would be produced in the longitude distribution. 
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Mills (1959) fits the Sydney 85 Mc/s observations of the continuum emission 
to a spiral arm model of the Galaxy, in which the spirals obey an equation of the 
form 


° 
R=0:'22R, exp |outec : 


7 


The parameters being chosen to give the best fit to marked ‘‘steps’’ in 
the longitude distribution of the 85 Mc/s radiation In the southern sky 
(270° <I! < 360°) the steps are well marked, and are such as might be expected 
in a spiral model of the Galaxy in which the magnetic field deviated randomly by 
about 20° from perfect alignment with the spiral arms. In the northern sky 
(0° </<go°) these steps are less well marked. 

The longitude distribution of radio emission in the northern sky is clearly 
shown in Fig. 1, in which profiles at constant latitudes are drawn for the three 
surveys used in the analysis. A critical examination of these figures shows that 
the distribution is complex even when contributions from easily recognized point 
sources have been removed. Since the profiles close to the galactic plane include 
both the thermal and non-thermal components, it might be thought that the 
profiles at constant latitude, one or two degrees from the galactic plane, would 
show the structure of the non-thermal radiation more clearly. Inspection of 
these profiles at + 2° and + 4° latitude (also Fig. 1) indicates no simple pattern in 
the emission. A striking step pattern is observed in the 408 Mc/s and 1390 Mc/s 
profiles at b1!= —2° (with steps at /7=10°, 30°). These steps do not agree 
well with the steps found by Mills (at /'= 14°, 24°, 40°), nor are they confirmed 
by the profile of b'=+2°. We conclude therefore that an examination of 
profiles of the northern Milky Way does not lead to a clear recognition of ‘‘steps’”’ 
in the longitude distribution, and that most of the features might be due to 
irregularities in the emission. 

The separation of the galactic emission into thermal and non-thermal com- 
ponents is illustrated in Figs. 2 and 3, and it is natural to look for evidence of a 
spiral structure in these diagrams. The isophotes of non-thermal emission in 
Fig. 3 are drawn from the profiles of Fig. 2 with, inevitably, a certain ‘‘ smoothing”’ 
of the data. We believe that all the significant large scale features of the radiation 
are preserved in this isophote diagram. We draw attention to two features of the 
longitude distribution of the non-thermal radiation. This first is a marked 
discontinuity at approximately /1=14° (the position of one of Mills’s steps). 
The second feature is a region in which the brightness temperature remains 
constant with increasing longitude between /1=20° and /#=40°, and then 
decreases rapidly, reaching a low value at /!!= 60°, the limit of the present survey. 
Beyond this limit and as far as /1=go°, the non-thermal radiation from the 
Galaxy seems to be weak, although a full survey at 408 Mc/s has not been made. 
The radiation from the Cygnus X region has been shown to be almost entirely 
thermal (Mathewson, Large, and Haslam 1960). 

These two features of the non-thermal emission might be fitted to a spiral 
model of the Galaxy of the kind discussed by Hanbury Brown and Hazard if ‘n’ 
is taken to lie between 1 and 2 and the direction of the magnetic field does not 
deviate randomly from the mean direction by angles which greatly exceed 20°. 
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There would be two spiral (or circular) arms of radii 2kpe and 6kpc, the line- 
of-sight being tangential to the arms at /T=15° and /4=50°. The positions of 
these arms are significantly different from those fitted to the thermal emission. 

6. Conclusions.—It will be clear from the foregoing remarks that any attempt 
to fit the continuum observations to a spiral model of the Galaxy is at present 
speculative. Our understanding of the emission will be much improved when 
the thermal (i.e. high frequency) component has been surveyed from the southern 
hemisphere with reasonably good resolution. Until such a survey has been made 
and until we have definite information on the position and direction of the galactic 
magnetic field we believe that identification of steps in the galactic radio emission 
with spiral arm structure is unlikely to prove correct. 

The present work has drawn attention to certain large scale features of the 
Galaxy seen from the northern hemisphere. These are the two thermal maxima 
observed at /¥=25° and /!!=80°, and the broad ‘“‘steps’’ in the non-thermal 
emission at /4=15° and /4=50°. ‘These features are well marked and can be 
interpreted in terms of a spiral model of the Galaxy, only if one postulates an 
independent distribution for the thermal and non-thermal emission. 
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ON THE COUNTING OF RADIO SOURCES 
IN THE STEADY-STATE COSMOLOGY 


F. Hoyle and }. V. Narlikar 


(Received 1961 June 13) 


Summary 


The problem of the number count of radio sources as a function of the 
incident flux is shown to depend on two crucial features: (i) the size and 
behaviour of condensations (ii) the dependence on age of the probability of a 
galaxy being a radio source. Provided the probability rises by a factor ~ 10° 
for galaxies with ages from H~- to about 2-5H-, and provided the primary 
condensations of the steady-state theory possess initial dimensions of order 
30 megaparsecs, the radio source count can rise more steeply than is the case 
for sources uniformly distributed in Euclidean space. Each primary con- 
densation contains of the order of 10° galaxies, which in the main expand 
apart from each other as the universe expands. 

It is shown that a luminosity function can be chosen for the radio sources 
giving results consistent with observation, not only for the source count, but 
also with the data on angular diameters and with experience in the problem 
of optical identifications. 





1. The log N—log P curve.—Six years ago, Ryle announced in his Halley 
lecture (16) that counts of the number N(P) of radio sources brighter than P did 
not behave in the manner to be expected of uniformly distributed objects in 
Euclidean space. 

From the point of view of the theoretician, the situation over the last six years 
has seemed uncertain, since different observers have been in disagreement over 
the real value of the slope of the log N—logP curve. In a comparison over a 
common portion of the sky, the slope obtained at Cambridge (1) was —2°5 at 
81-5 Mc/s and — 2-7 at 159 Mc/s, whereas Mills, Slee and Hill (2) obtained — 1-8. 
Very recently, the Cambridge observers have announced the result of a new 
survey, yielding about — 1-8 for the slope, and this they regard as superseding 
their former estimates. Another recent survey, at 960 Mc/s, by Kellermann 
and Harris (3) gives — 2-0. 

Although the concordance between the recent Cambridge work and that 
of MSH suggests that the slope probably is steeper than the Euclidean value 
of — 1-5 (a view expressed by Bolton in his survey given at the U.R.S.I. General 
Assembly, London, 1960 September), it must be noted that Mills has remarked 
on the likely existence of a systematic progressive error in his judgment of the 
values of P. In a ‘‘noise-limited’’ instrument, such as the Mills Cross, the 
observer tends to overestimate P more and more as the limit of detectability is 
reached. Mills has suggested that this effect could reduce the slope to — 1-65, 
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or even to —1°5 (4)*. On the other hand, the Cambridge instrument is not 
noise-limited, and no such systematic effect should be present. The concordance 
between the Sydney and the Cambridge work may therefore be coincidental. 

In the present paper we propose to accept the Cambridge claim that the slope 
is steeper than —1°5 for values of P less than 3 x 10-*° watt m~*(c/s)-!. We do 
not propose, however, to attach importance to the precise value of —1-°8. In 
fact, our best result for the slope will turn out close to — 1-6. We do not regard 
the difference between this value and the recent Cambridge result as a significant 
discrepancy. 

At P>3x10-**w.m~*(c/s)~! the slope may well be less steep. Indeed, 
Mills (4) gives a slope no greater than — 1-3 at high flux values. Certainly it 
does not seem necessary for the theory to provide a slope steeper than the 
Euclidean value at large P. (The values of P mentioned above refer to a 
frequency of 160 Mc/s. This frequency will be taken as standard throughout 
the paper.) 

2. The general problem.—Before the steady-state cosmology was proposed, 
ambiguities made the predictions of theory uncertain: three choices were 
available for the world curvature, and a double infinity of parametric values 
was also available for adjustment (the parameters being the cosmical constant 
and a constant of integration). Steady-state cosmology does not suffer from these 
ambiguities. There is no uncertainty concerning the line element, except for 
the numerical value of Hubble’s constant. Once this has been fixed from 
observation, purely geometrical questions can be settled explicitly. 

A further feature adds a still greater degree of definiteness. Because the time 
axis is open, into the past as well as the future, and because the line element is 
stationary (but not of course static), all observable properties must reproduce 
themselves, otherwise the continuing expansion would weaken and ultimately 
destroy the property in question. For example, galaxies must continue to 
condense, otherwise the spatial density of galaxies would tend to zero. Dynamo 
action must preserve any intergalactic magnetic field that may exist; and so on. 

These advantages seem to have encouraged an erroneous optimism that the 
theory is capable of handling explicitly a physical property whose nature is not 
specified—for example, that the theory will predict the slope of the log N—log P 
curve for radio sources without it being necessary to specify in any precise way 
the properties of radio sources. In a recent criticism of the steady-state theory 
Ryle and his colleagues (15) obtained a log N—log P curve without apparently 
needing to use any of the physical properties of radio sources. Manifestly, 
there is a paradox in this, for one cannot count objects without some close 
specification of what are the objects to be counted. The paradox is actually 
resolved by the circumstance that the mathematical analysis necessary to arrive 
at the log N — log P curve used by the Cambridge radio-astronomers itself imposes 
physical properties on the radio sources. If it should turn out that the properties 
thus imposed are indeed satisfied by the radio sources, then the claim that 
observation contradicts the steady-state theory may well be correct. 

In the following sections we shall construct a counter-example, showing that 
properties can be imposed on the radio sources that lead to a log N—log P in 


* We are indebted to Drs Allen, Palmer and Rowson for the information that the effect 
reported by Mills has been investigated also at Jodrell Bank. The effect is found to be present, 
and of quantitatively the order stated by Mills. 
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satisfactory agreement with the observations. It will be necessary to use a 
definite picture of the mode of formation of galaxies. Before beginning the 
cosmological discussion it will be useful to review briefly the astrophysical 
features of this picture. It is perhaps relevant to add that the essential features 
of the following section are not new (cf. (5, 6)). We shall therefore be using 
an astrophysical picture that arose quite independently of the present radio 
source problem, not a picture invented specifically for explaining the radio 
observations. 

3. Condensations in the steady-state theory.—It has been realized for some 
years that the formation of galaxies, and of clusters of galaxies, implies the 
existence of a distance scale apparently quite independent of the cosmological 
distance scale. In the steady-state cosmology, for example, the constants 
appearing in the cosmological equations are c, H, and G; and from these cH— 
is a length. But this is the so-called ‘‘radius of the observable Universe’’, 
about 3000 mpc (H=100kms~'(mpc)-'), very much larger than ~3 mpc 
given by the galaxies and their clusters. 

To obtain a distance scale appropriate to the galaxies and their clusters, 
it was suggested (5) that condensation may occur from an intergalactic gas 
at high temperature. Write V for the velocity of sound in the gas, then VH-! 
is a length. To obtain ~3 mpc, with H=100kms—(mpc)-!, V must be 
~300kmsec~!. It appeared immediately significant that this latter value is 
just of the order of the internal motions found in the galaxies, and also of the 
order of their peculiar motions. 

It was further suggested (5) that condensation is promoted by cooling in 
the hot gas—locally-cooled regions being compressed by surrounding hotter gas. 
The idea that pressure gradients are necessary for producing condensation has 
been supported by an investigation of Harwit (7), who has shown that in the 
steady-state model condensation cannot take place through the operation of 
purely gravitational forces. This is clearly connected with the remarks of the 
previous paragraph—if condensation through gravitational forces could take 
place, a distance scale for the condensation process involving only cosmological 
constants would be expected to turn up in the theory, and this we know to be 
impossible. A similar argument applies to other systems of cosmology; since 
the appropriate scale for the condensations cannot be built from the cosmological 
constants c, G, A, it is improbable that condensation can take place under purely 
gravitational forces. This argument is strongly supported by an investigation of 
Lifshitz (8). In other systems of cosmology it is difficult, however, to ascribe a 
high kinetic temperature to the intergalactic medium at the onset of condensation 
of the galaxies (a temperature of 10° °K is required at the expansion stage where 
the density falls below 10-*> gcm~*). 

Investigation of radiation by free-free electron-proton collisions showed that 
cooling can take place in a time scale of the required order, viz. the reproduction 
time of the steady-state theory, }H~-1, provided the gas density is not lower 
than ~10-*’gcm~*. Initially, the cooling is slow, the energy of contraction 
being mostly dissipated by radiation. A large measure of dissipation implies 
that a bound, stable structure will very likely be formed. However, after 
contraction to about 4 of the initial dimension—i.e. to about 1 mpc, the density 
rising meanwhile to ~ 10-**gcm~%—the cooling accelerates rapidly, and there- 
after only a comparatively small fraction of the energy of contraction is dissipated 
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as radiation. The energy of contraction then appears as dynamical energy, 
implying a fragmentation into sub-units, the sub-units being individual galaxies. 

Cooling within an intergalactic medium of temperature ~ 10’ °K and density 
10-*’ gcm~* therefore fitted the broad requirements of a theory of the formation 
of galaxies and of clusters. It gave the dimensions, time-scales, and fragmentation 
properties correctly. Moreover, rapid cooling by hydrogen ceases at about 
10°°K. The requirement that the hydrostatic pressure within and without a 
highly cooled zone be essentially the same, thus yielded a density ~ 10-** gcm-3 
inside a highly cooled zone, and this is precisely the order of the densities that 
must have occurred in the galaxies before flattening to a disk structure took place 
in many cases. If, further, we ask: what mass of gas, at density ~ 10-4 gcm~-%, 
can control dynamical motions of order 200-300 kmsec™, the answer is about 
10'Mo©. This was an additional cogent result—the theory immediately 
determined the general order of the masses of the main class of galaxy (smaller 
galaxies were thought to be splinters from the dynamical disruption of lesser 
masses). Finally, the hydrostatic pressure in question, ~ 107!*dynecm~, is a 
value that still affects the structure of our own galaxy in a remarkable degree. 

So far, the former theory appeared to be on a sound basis—the results obtained 
were sufficiently numerous to be reckoned a satisfactory return for the single 
important hypothesis of a high kinetic temperature (the assumption of a gas 
density ~10~*’gcm~* cannot be regarded as a serious assumption, since this 
is just the general order of the mean densities actually found observationally 
to be present in many clusters of galaxies). ‘To bring these considerations into 
closer relation with the steady-state cosmology it was found necessary to introduce 
further ideas, however, and these were more removed from the astronomical and 
astrophysical evidence. Accordingly, it is of great interest to find these further 
ideas turning out to be of particular relevance to the subject of the present paper. 
It does not seem likely that without these ideas the radio source counts could be 
explained within the framework of the steady-state theory. 

The mean intergalactic density appearing in the steady-state theory is of order 
10-*° to 3 x 10-8 gcm~%, the exact value depending on H and on the particular 
mathematical formulation of the theory. Hence 10-*’gcm~-% cannot be used 
everywhere in space but only in regions that are already partially condensed. ‘The 
question thus arose as to how such regions became partially condensed. 

It was suggested that condensation occurs in two stages, a primary stage from 
10-%gcm-* to 10-?7gcm~-%, and a secondary stage from ~10~-**gcm~-* to 
~10~-*4gcem~-%, the secondary stage being controlled by radiative cooling in 
the manner described above. Pressure gradients also play an important role 
in promoting primary condensation, the basic concept being one of equal 
hydrostatic pressure, the separate stages following the lines of the following 
table: 


Density Temperature Distance Mass 
(g cm-*) (°K) scale (mpc) (Mo) 
Before primary 
condensation as 10° 30 1016 
After primary 
condensation io 107 3 to" 
After secondary 


condensation 10-4 10! orl tol! 
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The differences of mass in the last column take account of the effects of 
fragmentation during condensation. Since the distance scales are not accurate 
within a factor ~2, the masses should not be taken as accurate within ~ 10, 
except in the last line of the table, where the mass is determined from the criterion 
of stability given above, not simply from multiplying volume by density. 

A primary temperature of ~ 10° °K implies an energy of ~ 2 x 10!" erg per gram 
of hydrogen. It was noticed that this is just the amount of energy to be expected 
if the newly created matter of the steady-state theory consists of neutrons, 
allowance being made for the energy lost by neutrinos in the neutron decay. 
The existence of a high kinetic temperature could thus be explained in terms 
of the form of the newly created material. 

The main question that now arose was: how does primary condensation take 
place? Not by radiative cooling; it was shown (5) that even under rather 
favourable conditions not more than about 5 per cent of the thermal energy can 
be dissipated by ordinary radiative processes. A possibility considered by Gold 
and Hoyle (g) was that, since the intergalactic material on the present picture 
possesses widely different temperatures in different places, ‘‘ heat engines ’’ may 
operate without the necessary dissipation needing to be even as high as 5 per cent. 
Such an engine could convert thermal energy almost entirely into the dynamical 
energy of mass motions of the gas. The precise nature of the heat engine was 
not specified, however, by Gold and Hoyle. It was thought that where the 
thermodynamic possibility exists for the operation of heat engines of great 
efficiency such engines are almost certain to arise. The situation has a logical 
analogy to the case of convection, where one considers the phenomenon to exist 
even though the detailed properties of the convective cells are not precisely 
specified. 

If the fraction of energy dissipated in primary condensation is small, thermal 
energy being converted largely into the dynamical motions of a number of 
condensed units (which subsequently undergo secondary condensation) then 
the units are most unlikely to stay bound together as an aggregate. Rather do we 
expect primary condensation to result in a number of clusters of galaxies in motion 
with respect to each other at speeds of the order of the speed of sound in a gas of 
kinetic temperature 10°°K, viz. ~3000kmsec~!. That is to say, we expect 
primary condensation to result in a number of comparatively dense blobs of 
gas (~10-*’gcm~%) moving with speeds ~3000kmsec™ with respect to each 
other. The largest of these blobs possess radii ~ 3 mpc, while the whole group 
of blobs possesses a radius ~30 mpc. For each blob of radius ~ 3 mpc there 
may well be many smaller blobs. 

Although the blobs possess kinetic energy large enough to expand the group, 
there is no reason why the initial motions should all be directed radially outward. 
Inward motions can occur, and may well be more common than outward ones, 
particularly if the system has experienced compression from the surrounding 
hot gas. But inward motions can serve only to delay, and not to prevent, the 
eventual expansion of the group—inward motions eventually become outward 
motions. ‘The situation is similar to an element of hot gas without bounding 
walls. Even without collisions between the constituent particles, the element 
soon begins to expand, and eventually all the particles move outward. In a 
similar way, the blobs of our group expand apart from each other, assisted in the 
cosmological case by the expansion of the Universe. 
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We may sum up the features of the picture as follows: 

(i) A primary condensation consists of an appreciable number of blobs of 
gas with density ~ 10-7’ gcm-%. 

(ii) All galaxies forming within the blobs of a particular primary condensation 
are age-correlated. 

(iii) The blobs move initially with criss-crossing motions, their velocities 
with respect to each other being ~ 3000km sec". 

(iv) After the initial criss-crossing, the blobs expand apart from each other. 

(v) The whole age-correlated region has an initial radius ~30 mpc, and 
may contain from ~ 10° galaxies. 

(vi) The largest of the blobs may contain upwards of 10% galaxies, but small 
blobs are probably much more common than large ones. 

The need for criss-crossing motions, followed by expansion, was noted in (5), 
where it was pointed out that an intergalactic magnetic field cannot be maintained 
unless some such phenomenon exists. 

We have now reached the stage where it is necessary to consider the relation 
of primary condensations one to another. This will be conveniently treated in 
the following section. It remains in this section only to point out that in special 
cases blobs with a primary condensation may possess small relative velocities, 
and may come close enough to enable them to resist the general expansion of 
the Universe. It is of interest in this connection that both Shane and Wirtanen 
(10) and Abell (11) find cases where clusters of galaxies form groups of 
dimensions ~20 mpc, a value satisfactorily close to the radii of our primary 
groups. 

(Although Abell gives a value ~ 20 mpc, this refers to Hz 180 kms~ (mpc)-}. 
For the distance scale of the present paper, H= 100kms~!(mpc)-}, Abell’s value 
should be increased to ~ 40 mpc.) 

4. The discrete model.—The existence of a condensation distance scale 
introduces a discrete element into the steady-state theory. In accordance with 
the discussion of the previous section, this distance scale will be taken as that 
associated with the primary condensations. 

To preserve the steady-state characteristic it is necessary that every feature 
of the Universe shall repeat itself in a time interval ~}H-1. The distribution 
of primary condensations existing at any moment is maintained, in spite of their 
expansion apart from each other, by the formation of new condensations. The 
interval }H~-! may be spoken of as the “‘length of a generation’’. 

The discreteness referred to above is inherent in the theory. For the sake 
of simplicity, we shall introduce a second discreteness representing a convenient 
approximation, recognizing however that this second discreteness would not enter 
intoastrict theory. We shall regard the inherently discrete primary condensations 
as arising at particular moments of time, the moments being separated by interval 
4H-1. Suppose that time ¢ is one of these moments. Then at time ¢ a system 
of primary condensations comes discretely into being. Between ¢ and t+4H-} 
no further primary condensations arise, the system of condensations born at 
time ¢ simply expand apart from each other in accordance with the de Sitter 
expansion factor R(t)=exp (Ht). But at time t+4H~— a new system of primary 
condensations again comes into being, the individual condensations and their 
mean spacing apart being exactly the same as was the case for the system born 
at time ¢. And so on for subsequent generations. 











Ree 














Aer UMN EIA HON 


[OU ET RIRRNT. oS 


BER e SERENA aS. Oe. ie nT Ee 








No. 2, 1961 On the counting of radio sources 139 





And we shall introduce the further artificial discreteness that the lattice formed 
by the primary condensations of each generation always has the same structure. 
For simplicity, we take the lattice to be cubic, and we also take the primary 
condensations to be spheres having centres at the lattice points of the cubic 
structure. After an initial delay we expect each primary condensation to join in 
the general expansion—once the delay caused by the criss-crossing motions 
discussed in the previous section is over, expansion of the lattice units keeps 
step with the expansion of the lattice itself. ‘That is to say, the radius of a lattice 
unit maintains a constant ratio to the lattice distance. We take the ratio to be }. 
This requires the lattice units of each generation to be well separated, since the 
choice of this ratio requires the distance separating two neighbouring condensa- 
tions A and B to be made up in the following way: 4 in A, 4 in neither A nor B, 
4 in B. 

Our choice of this particular ratio has been guided by the consideration that 
if the lattice units were initially almost touching each other, and if a time ~}H- 
be allowed for the delay in the expansion of the lattice units, then the ratio must 
be }e—3, which is close to }. 

Write a, for the radius of each lattice unit and /, for the lattice spacing at the 
beginning of uniform expansion. Then, from what has been said, a, ~ 30 mpc, 
l,~g0 mpc. The precise numerical values of a,, /, will not affect the following 
work in any important respect—their values could be changed by a factor of at 
least 2 without appreciable consequences. On the other hand, the ratio a,//, 
does turn out to be important. A decrease in this ratio would make the number 
count problem simpler, although a marked decrease would probably raise a 
different difficulty in that non-isotropic effects would then arise in marked degree. 
The latter question will be taken up later in the paper. 

It is convenient to describe the lattice characterized by (a,, /,) as the “‘ first 
generation ’”’ lattice. This will be the ‘‘ youngest ”’ lattice that an observer will 
expect to find. The observer will also detect older lattices. ‘Those older by 
(n—1)/3H, n=2, 3, ... will possess lattice units and lattice spacing characterized 
by (a,,, /,,) given by 


(4,,, f,) =e"*8(ay, 4). (t) 


We shall refer to (a,,/,) as the ‘‘nth generation lattice’’. Allowing for the 
initial delay period of 4H-1, the age of the mth generation is just 4nH-. 

The galaxies to be found in the lattice units of a particular generation will not 
possess ages as long as that of the lattice in question, since an interval of time 
must elapse before the galaxies and the stars within them are formed. We 
tentatively set this interval as ~#H-1. Then galaxies containing stars appear 
in the second generation. Galaxies possessing stars of age ~ H—! appear, however, 
only in the 5th generation. The 5th generation lattice has a,;~ 100 mpc, 
l;~ 300 mpc, although once again, we note that these numerical values could 
easily be in error by a factor ~2. 

Since the majority of definitive identifications of radio sources are with 
giant E galaxies, and since the stars of E galaxies give a late colour index, and are 
almost certainly at least of age H-1 (12), it appears that radio sources occur 
mainly in old lattices, in generations later than the 5th. We therefore see that 
the discreteness inherent in steady-state cosmology must enter the radio source 
problem on a scale at least as great as /;~ 300 mpc. For the most part in what 
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follows we shall be able to employ simple averages derived from the usual 
continuous form of the steady-state theory. ‘The present discrete picture is 
not then needed. At certain stages of the argument, however, the discreteness 
of the condensation picture is required, indeed at one crucial point. It is important 
to notice that this crucial aspect is concerned with the inherent discreteness of the 
condensation picture, not with the artificial discreteness of the sequence of 
regular lattices. 

Several further points will be used at later stages. ‘The probability that an 
observer at random lies inside a lattice unit belonging to the mth generation is 


given by 
‘, (7) = “ (i) = as 2) 
j 3 In 3\4/ 64 


the ratio 4 being used for a,/l,. The probability is therefore the same for all 


generations. 
Although f is appreciably less than unity, the observer must lie in lattice units 
belonging to an infinite sub-set of nm=1, 2,.... But since the space density of 


galaxies weakens with increasing m as exp(—™m) the main contribution to the 
density of galaxies at the observer’s position comes from the smallest value of n 
in the sub-set. 





1sT 
4TH GENERATION. GENERATION 


Fic. 1 


The oldest stars in our Galaxy appear to have ages of order 
1-5 x 10° years~ 4H-1/3. This implies that we ourselves are inside a sixth 
generation lattice unit. We may also lie inside a younger lattice unit, and we 
must certainly lie inside older systems. The latter make no more than a small 
contribution to the total density of galaxies, while galaxies belonging to any 
younger lattice may still be largely undeveloped. In this connection we note 
that our presence here on the Earth guarantees some such situation, since our 
presence demands a galaxy, stars, planetary formation, etc. 

It has sometimes been objected to the steady-state theory that the mean age 
of all galaxies should be }H-", whereas in our locality observation suggests ages 
that are generally of order H-!. We now see that our presence inside a fifth or 
sixth generation lattice unit would give systematic ages ~ 1-5 x 10!° years for 
the developed galaxies around us. 











No. 2, 1961 On the counting of radio sources 141 


In Fig. 1 we show the forms of lattices for the first and the fourth generations, 
the scale in the latter case being greater by e. The positions of the centres of 
the lattice units of any particular generation are probably anticorrelated with those 
of preceding and succeeding generations. This is because the formation of a 
primary condensation at a particular place probably inhibits the formation of a 
condensation at the same place in the succeeding generation. Since fis only ~4 
there is opportunity for the lattice units of any generation to form in the interstices 
of the lattice of the preceding generation. 

5. The general mathematical formulation.—The line element of the steady-state 
theory is 

ds* = c* dt? — exp (2Ht)[dr? + r?(d@* + sin? 6 d¢?)], (3) 


where r=0 is the observer’s position. Our discussion of the lattice distances, 
given in the preceding section, refers to a particular value of t, distances being 
given by exp(Ht)j{du. The line integral with respect to u is here to be taken 
along a straight line in the Euclidean sub-space 


du? = dr* + r?(d6? + sin? 0 d¢?). (4) 


Indeed, if we choose the arbitrary zero of ¢ to be the present, then the present 
lattice distances are simply coordinate distances. 

The lattice points of any generation possess coordinates 7, 0, d independent 
of t. And once the lattice units reach the stage of uniform expansion, the clusters 
of galaxies that make up each lattice unit may also be taken as having constant 
coordinates, if we neglect the finite size of a cluster. 

Taking t=o as the present, we ask: what is the r-coordinate of a galaxy 
that we observe (at the present) to give a red-shift AA/A? The answer in the 
steady-state theory is 


r=cH—2, (5) 
where 
AA 
=_—, 6 
s=5 (6) 


A galaxy of intrinsic bolometric luminosity L and red-shift z will be observed as 
having an apparent bolometric intensity P given by 


LH? 


O° adr ta * 


where P is the energy flux normal to the direction 0, ¢ of the galaxy in question. 

Suppose now that we wish to count the number of galaxies brighter than some 
specified P, on the basis that all galaxies possess the same intrinsic L. First we 
define a value of z, z,,, say, by solving (7) for z, using the specified P. Our problem 
is to count the number of galaxies with z<z,,. We attempt this, first using the 
simple continuous steady-state theory. In this the number of galaxies per unit 
proper volume is a constant, independent of the epoch. We consider observations 
made at t=o. The proper volume in the coordinate range r to r+dr is 4nr* dr 
at t=o. The number of galaxies in this coordinate range is therefore 


(constant) x r? dr. 


By (5), the coordinate range r to r + dr corresponds to a red-shift range z to z + dz, 
where z is given by (5), and dz=c-!Hdr. ‘Thus the number of galaxies giving 
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red-shifts in the range z to z+ dz is of the form 
(constant) x 2? dz. 


Hence, at first sight, the answer to our problem is given by integrating the latter 
expression from z=o0 to z=2,,, viz. (constant) x2z,,3. Since the constant is 
independent of the epoch, and also of P, a similar calculation for varying P 
gives the number count as a function of P. The count is not absolute, of course, 
since we have not specified any absolute proper density of galaxies. 

This simple procedure overlooks, however, the point that some of the galaxies 
we have counted, and which now exist at t = 0, did not exist at the time the radiation 
would need to have been emitted in order to reach us at the present day. For 
red-shift z the radiation would need to have been emitted a time H-In(1 +2) 
ago. 

To cope with this complication we must introduce the age distribution of the 
galaxies. In the continuous steady-state theory the number of galaxies with 
present-day ages between 7 and 7+dr is 


(constant) x exp (— 3H7r) d7r per unit proper volume. 


Thus if we are considering galaxies giving red-shifts between z and z+dz we 
must restrict ourselves to those in the age distribution with 7>H-1In(1+2). 
This number is 
(constant) x 2° dz | exp (—3Hz) dr, (8) 
H-"in(1+2) 
where the upper limit of the integration with respect to 7 remains to be stated. 

The three-dimensional sphere at t=o defined by a red-shift z, or by the 
corresponding r=cH~-1z, must now be compared with the lattice spacings of 
our model discussed in the previous section. If for generation n the lattice 
distance /,, exceeds 2r then certainly not more than one lattice unit of the mth 
generation can fall completely inside the sphere. The phenomenon is then 
essentially discrete. Continuous averages are no longer applicable. 

The point of transition from the continuous picture to the discrete model 
is of necessity rather uncertain. We could, for instance, require that at least 
two adjacent lattice units of generation m should fall completely inside the sphere 
of radius 7 in order that the continuous picture be applicable up to the age implied 
by generation n. This requires /,, < 6r/5. 

To deal with the slight measure of uncertainty involved in the transition to 
discreteness we apply the continuous picture up to a generation m defined by 


l, =kcH—z, (9) 


where k is a constant between 1 and 2. Ina later section we shall show that the 
precise value of k is not important. 
It remains to relate /,, to the lattice distance of the first generation. We have 


=~ a—I 


1, =l, exp =cH-1z, exp 








With /, and z, specified for the first generation, the limit of applicability of the 
continuous theory occurs at generation n, determined by solving 


=f 


3 





(10) 


kz~2z,exp 
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The age of the mth generation is 4nH-1. By (10) this is 


H>(3 +1n =). 


However, this does not give the correct upper limit of the integral in (8), since 
the age of generation m has been reckoned from the time of formation of the 
primary condensations, whereas the zero of 7 in (8) refers to galaxies of zero age. 
There is an interval of time, A say, between the formation of a primary condensa- 
tion and the formation of galaxies within it. Above we took A to be ~§$H™. 
The correct upper limit is thus 


-(isinl2)—a, 


a 

which can be written as H—! ln kz/z,', where 
2,'=2,exp (HA—}4). (11) 
Thus the count for the range z to z+dz, given by the continuous picture, is 


kz 


-1 in & 
(constant) 2? de * exp (—3H7) dr. (12) 
H-* In(1+2) 
The transformation y= Hr—In(1+2) reduces (12) to 
(constant) . aol exp (— 3x) dx. (13) 


To obtain the total count down to flux level P, i.e. for z up to z,,, we integrate 
(13) as follows 


kz 
tm og. Oe 
7 | O+94" exp (— 3x) dx. (14) 


“x (1+) 
k—z,’ 


(constant) : | 


é 


The integral is taken over an area in the (y,z) plane enclosed by the z-axis, the 
line s=<,,, and the curve 
kz 
x= nea (15) 
This gives the number of galaxies brighter than P, in so far as this can be determined 
from the continuous theory. 

Over and above (14) there is the possibility that a further contribution arises 
from discrete effects. These do not have much importance, except in the case 
where a sphere of radius r, r=cH-1z, (z<z,,), is invaded by a lattice unit of a 
generation even older than the value of m corresponding to z=2,, in (10). Such 
a case represents a genuine fluctuation that must be given special consideration. 
We shall return to this question at a later stage. 

It is already a matter of interest that (15) vanishes for 





2,’ 2,’ d ; 
2, = age knits and 2,/<€1). 
m k—2,' k ( 5 1 < ) 
This implies that as z,, increases above ~2,'/k the product NP%? increases to 


a maximum value and thereafter declines. Hence we have a log N—log P curve 
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that has a slope initially greater than —1-5. This is not our answer to the 
criticism of Ryle, since the increased slope only applies for z,, somewhat greater 
than 2,'/k, which will turn out to be about 1/15, whereas the observed slope 
maintains its steepness up to z,,~0°3 or 0-4. Ryle has argued that the slope 
is maintained up to z,, = 1 or more, but it seems to us that the evidence concerning 
angular diameters of radio sources scarcely supports this extreme view. 

The simplest hypothesis capable of steepening the log N—logP curve at a 
larger value of z,, is the following: only galaxies of age greater than 7, say, 
become strong radio sources. We must then count the sub-class of galaxies 


with age greater than 7. The number of these galaxies varies with z,,, and hence 
with P according to 
kz 
‘2m g2dz In ———=7. 
(constant) ——,| “** exp (— 3x) dx. (16) 
(1+2)9J ar 


In order that the integral with respect to x shall yield a positive result it is necessary 
that 

——— _ >exp (HT). I 
Gta or ) (17) 
For this to be certainly the case the lower limit of the integral with respect to z 
must be taken at the value of z given by considering the equality sign in (17). 
The count can only be applied for values of z,, greater than this particular value 
of z. Indeed, the count falls to zero at z,, equal to this z, and the log N—log P 
curve is steeper than —1-5 for z,, somewhat greater. That is to say, the slope 
exceeds — 1-5 for z,, somewhat greater than the value of z given by 


kz 
(1+2)2,' 
For T between 2H~ and 3H- this yields a value of z of ~0-3. 

These considerations are a special case only, however, of the general situation 
that arises when the probability of a galaxy being a radio source is a function of 
the age of the galaxy. We consider the problem in the following terms. 

(1) It may be that only galaxies with some special property become radio 
sources of the intense type that contributes appreciably to the log N—logP 
curve. Examples of such properties might perhaps be those galaxies that are 
members of physically connected pairs, or perhaps E-type galaxies. So long as 
the property is not age-correlated, however, we simply choose out the sub-class 
of galaxies that possess the property, or properties, in question. This changes 
the counting scale but not the relative counts for variable P. 

(2) The intrinsic radio luminosity will certainly not be the same for all radio 
sources. So long as the intrinsic luminosity is not age-correlated, however, it 
is possible to represent the luminosity function by a histogram, considering each 
section of the histogram in the manner described above, each section leading to 
an expression of the form (14), but each with a different z,, (for the same P). 
Varying P gives a log N—log P curve for each section of the histogram. If none 
of these separate curves possesses slope greater than — 1-5, the combined curve 
cannot possess slope greater than —1°5. ‘The luminosity function can be 
represented to any required degree of accuracy by taking a sufficient number of 
sections in the histogram. 


=exp (HT). (18) 
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(3) The situation is changed only if either the probability of a galaxy being 
a radio source or the luminosity function is age-dependent (or of course both 
of these). ‘The dependence of the probability on age we denote by K(r). To deal 
with an age dependence of the luminosity function, we represent the luminosity 
function by a histogram at each value of 7, the histograms being chosen in the 
following related way. ‘The luminosity function is first normalized to unit area 
for each 7, and the histograms are then made up of the same number of sections 
of the same areas. At different 7 the sections of the histograms are related serially, 
the first to the first, the second to the second, and so on. In the following 
considerations we consider only one set of related sections of the histograms, and 
we denote the dependence of the intrinsic radio emission for this section by L(7). 
The results for other sets are similar to those obtained below. 

Consider galaxies with present ages between 7 and ++dr giving red-shifts 
between z and z+dz. The number of these galaxies that are radio sources is 
given by an expression of the form 


(constant) 2? dzK(7— H-!ln1+2) exp (—3H7r) dr (19) 


provided the values of z, 7 can be related within the terms of the continuous theory. 
The relevant condition is that already discussed above, that 


Hr<in & ; (20) 
%y 
Introducing y= H7—In(1+2) (19) and (20) can be written as 
# , 
(constant) Gta) K(x) exp (— 3x) dx dz, (19) 
kz 
< In ——_.. F 
i (1+2)2,' ails 


A second condition is imposed by the requirement that the apparent bolometric 
flux be greater than P, viz. 
H*L(x) 


~ 4mc?z?(1 +2)?" 


(21) 


The radio astronomer does not of course measure the bolometric flux, but 
the flux per unit frequency range at some specified frequency. However, it can 
readily be shown that for sources with frequency spectra dv/v the measurement 
of flux per unit frequency range yields the same relative counts of sources as 
does the total bolometric flux. Most radio sources possess spectra of about the 
form dv/v°"8. This we regard as being sufficiently close to dv/v for the difference 
to be neglected. The effect of dv/v°*, and of regarding P as the flux per unit 
frequency range, is to replace the factor (1 +2)? in (21) by (1+2)'*. We have 
not troubled to take account of this slight modification. 

To obtain the number of sources with flux greater than P, (19’) must be 
integrated over the portion of the y,z plane defined by the conditions (20’) 
and (21). The (x, 2) plane is shown in Fig. 2. First, draw the curve 


x =Inkz/(1+2)2,’. 


Next draw the curve along which 


2(1+2)= ap x) (22) 
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The relevant area is bounded by these two curves and by the z-axis, in the 
manner shown in Fig. 2. 

It will be noticed that the age dependence of L does not affect the weighting 
function (19’), although it does affect the area of integration. The two functions 
K(x), L(x) do not therefore affect the counts in the same way. Indeed, it turns 
out that while K(x) can have an extremely important effect on the log N—log P 
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curve, the dependence of the luminosity Z on age is much less effective. This 
is demonstrated by an explicit example given in the Appendix. Since the form 
of L(x) is comparatively unimportant, it will be sufficient in later work to take L 
independent of yx, keeping only the probability function K(x). In this case the 
general integral reduces to a form very similar to (14), viz. 


2m 22 dz ER A )d ( ) 
—- : ex a ° 2 

marl, (x) exp (— 3x) dx 3 

k-z,’ 


(constant) | 


Writing g=exp x gives the somewhat more convenient expression 


kz 
| KE. (24) 


™m 
a (1+2)5), 
—2,’ 


(constant) | 


i 

6. The basic log N—log P distribution.—In order that (24) give a result of an 
essentially new character, it is necessary that an appreciable contribution to the 
integral with respect to q shall come from the upper part of the range of g. This 
will be the case if K(q) varies with g according to some power >3. The simplest 
case is where Kcq*, when the integral is equally weighted over the whole range 
of g. We shall adopt this distribution throughout the following work, noting, 
however, that essentially similar results are obtained for any power law between 
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the | about 3 and 5. In this connection it may be noted that it is always possible, 
over a limited range of g, to represent any general function K(q) by a power law. 
Thus Kocg* simply means that over the range of q in question we regard K as 


tin : 
i being represented by this particular power law. It will turn out that the relevant 
rns range of g covers the generations of our lattice model from n=5 to n=10, and 





e P } that values of z from about o-1 to about 0-3 are mainly involved. ‘The corre- 
sponding range of q is from 1 to about 3. Our dependence Kocq* therefore 
implies that the probability of a galaxy (perhaps of a certain restricted type 
—e.g. members of pairs) being a radio source increases over the relevant 
generations by a factor ~ 3%. A fuller discussion of these questions will be given 
at a later stage when numerical results have been derived, and when the various 





} other physical factors of relevance have been discussed. 
With Kocq', (24) gives 
( ‘ m «= dz ke | 
(constant) i ee tees I (25) 
k—z,’ 
} An important point now emerges. We shall be concerned mainly with values of 


z,, in the range o-15 to 0-3. Since the integrand of (25) is heavily weighted 
towards its upper limit, we can consider the factor 
kz ; 
(1+2)z,’ 
for values of z near z,,. It is seen then that kz/(1+2)z,' appreciably exceeds 
unity (k we expect to lie between 1 and 2, and k/z,' will be found to be ~15). 
This means that, without severe approximation, we could neglect the —1 in 


this factor, in which case the uncertain factor k simply goes outside the integral 
and is absorbed into the multiplying constant. Hence we see that the choice of k 


"his does not affect the log N~logP curve in any really important degree. 

rm To obtain numerical values, k/z,' must be specified. The value of z,' is related 
eL to the first generation lattice by 

ad 3,’ =2, exp (HA- }), 


: where 2,=cH~l, ~ 0-03 for /,~ 90 mpc (H=100 kms~!(mpc)~") and A is the 
time delay between the formation of a primary condensation and the development 
23) of the first radio sources within it (more precisely of the first radio sources that 
make an appreciable contribution to the count). The optical identifications 
available show a marked association of radio sources with galaxies of late colour 
: index, the stars of which can scarcely have ages less than H~!. ‘Thus, allowing 

3H~-! for galaxy formation, together with a further H—! for the ageing of the stars 

24) we have z,'~2,exp $=0'11 for z,;=0°03. The value given to k defines the cut-off 
from the continuous theory. We have seen that a value between 1 and 2 would 
seem appropriate. With k~1-6 2,’~o-11, k/z;'~15. This will be the value 











, “a used in succeeding work. As already emphasized, the results are not sensitive 
bea to this particular choice. 
; en The integral (25) can be evaluated explicitly to give 
es 
_*%1)\) ~) 4, Be +28 8 1128+ 1527+ 6a] %m 6 
(: ?) saiailiid 2i1+z? k 6(1+2)* a (26) 
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Taking k/z,’=15 counts on an arbitrary scale can be worked out numerically as 
a function of z,,. Since Pocz,,-#(1+2,,)~* the corresponding energy fluxes can 
also be given on an arbitrary scale. (It will be recalled that so far we are only 
considering sources of a single fixed intrinsic luminosity.) Results are given in 
Table I. The associated value of NP*, again on an arbitrary scale, are also given 
in Table I. 
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FIG. 3 


The log N—log P curve is plotted in Fig. 3. It is evidently steeper than — 1°5 
over the range from z=o0°1 to z=0°3. This range of z corresponds to a change 
by a factor 125 in P, as compared to about a factor 15 for the results upon which 
Ryle based his criticism of the steady-state theory. The slope of the log N—log P 
curve is close to — 1°8 for z between 0-2 and 0-3, and the slope even exceeds —2 
for z between o-1 and 0-2. 

To relate these results to observation it is necessary that the hitherto arbitrary 
scales of N and P be fixed so that N and P are in reasonable relationship to the 
values obtained in the observational surveys. To give a suitable association we 
take the unit of P to be 10-2 w.m~*(c/s)-! and we regard N as the number of 
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sources per steradian. ‘This transfer to absolute scales yields information 
concerning absolute values of the probability function K and of the intrinsic 
emission L. We shall defer a discussion of the probability function until a later 
section. 

The value of ZL can now be worked out from (7), by inserting z=0°3, 
P=2x10-*%w.m-*(c/s)-!s-1. The result is L=3-3x10%w(c/s)-!, the 
values of P and L being taken at ~160 Mc/s. For a source with spectrum 
dv/v°* extending up to v=3000 Mc/s, the total energy emission would be 
~ 4x 10" logs~}, a value that falls within a factor of about 2 of the total emission 
of the source Hydra A. It is of course the case that a value must be ascribed to H 
in (7). In the present calculation H = 100 kms~!(mpc)~! was used. 

7. Fluctuations in the number count.—Two quite different types of fluctuation 
enter the present theory: fluctuations in the number count due to the presence 
of the observer in or near one or more lattice units of the discrete model, and 
fluctuations from isotropy in the distribution of sources on the sky. The two 
types of fluctuation are affected by the discrete properties of the model in markedly 
different ways. The effect on the count is more extreme when the observer lies 
inside a lattice unit. ‘This lifts the log N—logP curve at large values of P, 
thereby reducing the overall slope of the curve to a greater degree than if the 
observer were to lie immediately outside a lattice unit—in the latter case the curve 
is unaffected at both large and small values of P, but receives a wiggle at 
intermediate values. Departures from isotropy only become an issue, on the 
other hand, when the discrete lattice units are projected on to small areas of the 
sky. Hence departures from isotropy must be considered for the case where the 
observer is distant from a lattice unit, not when he is inside it. The question of 
non-isotropy will be considered later in the paper. For the present, we are 
concerned with effects on the number count, particularly with placing the 
observer inside a lattice unit, or units. 

First we note that an observer cannot lie in more than one lattice unit of a 
particular generation. If he lies in two lattice units of different generations, 
one of the units will have appreciably greater weight than the other. Thus if 
the two units belong to generations with ages great enough for them both to 
belong to the range of q in which K(q) is rapidly increasing, but not so great 
that K(q) has reached an ultimate saturation level, then the older of the two units 
contributes the greater weight. It follows that we need only consider the presence 
of the observer inside one lattice unit, the one giving the greatest weight. 

The lattice units are never so large that the 1+221 cannot be used as a 
tolerable approximation for the sources within the particular lattice unit in which 
the observer happens to be situated. These particular sources, taken by 
themselves, therefore follow a log N—logP curve with slope close to —1°5. 
The general effect is shown by Fig. 4, where log NP?" is plotted against log P. 
Curve I is simply a plot of the values of Table I. At large P, the curve falls 
steeply away. Curve II of Fig. 4 shows the effect of a fluctuation. At large P 
the fluctuation dominates and NP?’ is effectively constant. Then as P decreases 
the curve first falls to a minimum, and thereafter rises steeply. The minimum is 
caused by reaching the limit of the local lattice unit at log P~ 2-0, the contribution 
of the lattice unit ceasing at this value of P. 

The reason why the fluctuation has been set at the level shown in Fig. 4 
will now be explained. We have seen that the age of our own Galaxy suggests 
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that we are situated inside a lattice unit belonging to about the 6th generation. 
The 6th generation lattice has spacing /,exp 5/3. ‘The distance of an observer 
at the centre of a unit cube of this lattice from each of the nearest eight vertices is 
(4/3/2)l, exp 5/3. This distance corresponds to a red-shift (4/3/2)c—H1, exp 5/3. 
With 1, ~ 90 mpc, H = 100km sec! (mpc) -1 (we have used these values through- 
out), the resulting value is 0-137. According to Table I the value of N at z= 0-137 
is ~17. If the whole of this count were contributed by the 6th generation alone 
then each of the eight lattice units of the unit cube would contribute ~ 2 to the 
count. This would ignore, however, the contribution of the 5th generation, 


so 
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and also the contribution from one or two lattice units of the 7th generation 
that may be observable at red-shifts <0-137. We have made allowance for this 
by taking half the count as coming from the 6th generation, in which case each 
of the eight lattice units of this generation contributes ~1 to N. 

The radius of a 6th generation lattice unit ~ }/,exp 5/3. A source at such a 
distance has red-shift ~ 4z, exp 5/3 (,=0°03). Thus an observer situated at the 
centre of a 6th generation lattice unit will find that all the local sources coming 
from the unit in question give red-shifts z <o0-o1 exp 5/3 =0°053. The value of 
log P at z=0°053, on the scale of Table I, is about 2-0. ‘Thus with N~=1 at 
log P = 2:0, the constant level of log NP? to be used for curve II of Fig. 4 is ~ 3:0. 

The precise depth of the minimum of curve II depends on the steepness in 
the region log P~1-5 of the fluctuation-free curve. This is affected by our 
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choices for z,, k, A. Since our choices are somewhat uncertain, the depth of 
the minimum, and its position with respect to P, is not well determined. All we 
can properly assert is that a minimum is to be expected in the region from 
log P= 1°5 to log P~ 2:0. 

If 10-6 w.m~* (c/s)? is the unit of P, then the values of N in Table | at 
P~1o-*w.m~*(c/s)~! are approximately correct for the number of radio 
sources counted per steradian. On this basis N~ 1 implies only about one source 
per steradian. It is clear, therefore, that although the fluctuation appears 
impressive when plotted in the form of Fig. 4, the actual number ot sources 
introduced by the fluctuation is very small compared with the count at 
P<10-*w.m~*(c/s)-+.. The precise height of the fluctuation of curve II is 
not therefore a matter of great importance. Doubling the height (i.e. increasing 
log NP*!? by 0-3) would only add about 1 source per steradian at flux levels 
>10-*4 w.m~*(c/s)~}, and no profound attention would be paid to this. The 
existence of the minimum is of considerable theoretical interest, however, since 
it clears up an apparently puzzling point. 

Consider a single object, with any probability function for being a radio 
source at rest in the usual cosmological frame of reference (r, 6, 6 coordinates 
constant). Suppose that observations of the source are made by an ensemble 
of observers at rest with constant r, 6,4 coordinates. Imagine the observers 
collate their experiences and that they draw an average curve relating the 
probability of observing the radio source against flux P, the averaging being 
performed with equal weight per unit proper volume for observations made at 
the same value of ¢. Then the curve so obtained will be the same as the 
log N — log P curve determined by a single observer in the case where an ensemble 
of objects, each exactly similar to the original object, is uniformly distributed 
throughout space. Hence an arbitrary set of objects distributed in space in any 
fashion will yield exactly the same result when averaged w'th respect to an 
ensemble of observers as will a grand ensemble of objects uniformly distributed 
in space. Now the continuous steady-state theory treats just the case of a grand 
ensemble—it treats a mixture of ages in the correct proportions. Taking account 
of the eventual saturation of any probability functions (so that divergences 
cannot occur) the smooth continuous theory yields the form of log N—log P 
curve used by Ryle and his colleagues. Hence the latter curve must represent 
the average arrived at by the ensemble of observers. ‘The discrete model of the 
present paper can lead to a particular observer obtaining a different log N — log P 
curve, but it cannot prevent the average with respect to all observers from being 
of the usual form. ‘The averaging of the observers must take place through a 
volume with dimensions greater than the lattice spacing of the oldest generation 
that makes an important contribution to the count. This is determined by the 
saturation of the probability function. The results derived above are unaffected 
by a saturation at the roth generation. The relevant lattice spacing could 
therefore be 1,exp3~20/,=1800mpc. Averaging must therefore take place 
over the whole ‘‘observable’’ region of the Universe (this is defined as 
cH-1~3000mpc, for H=100kms~!(mpc)-'). While it would be possible to 
raise logical objections to observers actually collating their experiences over 
regions as large as this, one can still carry through a mathematical average. The 
question arises as to how the results derived above could lead, when averaged 
for all observers, to the log N—leg P curve given by the usual continuous theory. 
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The answer lies in the minimum shown in curve II of Fig. 4. Some observers 
will be situated inside late generation lattice units, e.g. the 8th and gth. The 
fluctuation level of NP*? for an observer in such lattice units is not only lifted 
but the fluctuation extends to smaller values of P. The minimum of curve IT 
of Fig. 4 is shifted left by ~1-oin log P. ‘The observer would find a log N—log P 
curve with slope —1°5 at high P, the slope would become somewhat less than 
—1°5 as P decreased (due to the usual red-shift effect), but then for P around 
10-*° w.m~* (c/s)~+ the slope would fall precipitately away to the minimum—the 
slope would be markedly less than that given by the completely continuous theory. 
Eventually at low values of P the slope might increase somewhat, but NP*? 
would not rise back to its value at high P. 

We expect observers of the sort described in the previous paragraph to have 
a large weight in determining the average for the whole ensemble of observers. 
This is because such observers add comparatively large numbers of sources at 
high values of P. ‘The slope of a log N—log P curve can be appreciably changed 
by adding only a moderate number of sources provided these are added at high 
and at medium P values. No great differences in the counts at low P values are 
needed. Indeed we expect all observers to obtain much the same counts at low 
flux values. 

Can any explanation be given of why we do not ourselves lie inside a late 
generation lattice unit? On the basis of sheer probability the chance against 
lying inside a lattice unit of any particular generation is about 5/6 (cf. Section 4). 
The chance against lying inside lattice units belonging to the 7th, 8th and gth 
generation would thus be ~ (5/6)%, and this is not small. Moreover there could 
be biological and physical reasons of relevance, relating to our presence on the 
Earth. For life to be possible (in our particular form) rather stringent chemical 
conditions, as well as astronomical conditions, must very likely be satisfied. 
It is entirely possible that these conditions are only met with in galaxies of about 
the age $H~', in which case our presence inside a 6th generation lattice unit 
would not be an accident. And this could inhibit our chance of lying inside a 7th, 
or even an 8th, generation unit, since the lattice positions of succeeding 
generations are very likely anticorrelated. Nor is it an entirely unreasonable 
speculation that the degree of electromagnetic activity associated with 8th and 
gth generation lattice units might be so great that biological development would 
be significantly affected by a very high cosmic ray background inside such units. 

Quite apart from these points, however, it is important to notice that although 
a discrete model does not change the average of an ensemble of observers, 
it does make each observation yield a log N—log P curve that is different from the 
average curve. ‘The situation is not similar to the case of a Maxwellian distribution 
of particles, where most particles have energies close to the mean energy. A typical 
observer does not obtain a log N—log P curve close to the average of the whole 
ensemble. Observers obtain widely differing curves depending on their 
individual relationships to the discrete structure of the condensations. Only 
after a very large-scale averaging process has been carried out does the mean 
log N —log P curve appear. 

Our answer to Ryle’s criticism lies at exactly this point. Every observer 
taken individually obtains a log N—logP curve differing from the mean curve. 
Hence every observer could interpret his result as disproving the steady-state 
theory, even though in the mean their results actually verified the theory. For 
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example, an observer situated inside a late generation lattice unit (8 or g) could 
argue that as P decreased below 10-*° w.m~? (c/s)~! the slope differed markedly 
from the continuous steady-state curve, but in just the opposite sense to that of 
the terrestrial observations. We are dealing with a property in which the 
fluctuations make an important contribution, particularly fluctuations arising 
from old lattice units. The property is due to the weighting effect of the probability 
function K(7) which causes the mean situation to have no close resemblance to 
the situations in particular cases. Such an effect does not arise when K is 
independent of 7, as it is in the simple counting of galaxies. 

8. The effect of the luminosity function.—In the previous sections we have 
shown how a log N—log P curve appreciably steeper than — 1-5 can be obtained. 
In the present section we shall be concerned less with the precise value of the 
steepness of the log N—logP curve than with the discussion of a realistic 
luminosity function, and with the relation of the luminosity function to an absolute 
scale for P. To obtain the best fit of the theory to known data concerning the 
luminosity function some sacrifice of steepness has to be made—the value obtained 
below for the slope resulting from the combination of a number of luminosity 
classes will turn out close to — 1-6. 

We have seen in Section 6 that the unit of the logP scale in Table I and 
in Figs. 2, 3, 4 can be taken as 10-*°w.m~*(c/s)~! if the intrinsic luminosity 
of the sources in question is ~3 x 10%6w.(c/s)-! at frequency 160 Mcs-}. 
Such sources we refer to as of class I]. We consider five other classes: 
I, L=10% w. (c/s)-?;_ III, L=107' w. (c/s)-?; IV, L=3 x 107 w. (c/s); 
V, L=10*% w.(c/s)-!; VI, L=3~x 1078 w. (c/s)-1. Each class of source 
is taken as possessing a log (NP?’)—logP curve of the form of curve II of 
Fig. 4, but with an appropriate displacement of the flux scale. Thus class | 
is displaced by ~o-5 to the left (with respect to Fig. 4 for class II), class III 
by ~o:5 to the right, class IV by 1-0 to the right, class V by ~1°5 to the right, 
and class VI by 2:0 to the right. The maxima for I and II were taken as of 
equal height in the log NP?—logP plane, the maximum of III was reduced 
by o-2 in the logarithm, that of IV by 0-44, that of V by 0-24, that of VI by 0-64. 
‘These choices were made in order to give what we feel to be reasonable agreement 
with all the available data, such as it is known to us. The result of combining the 
six classes is shown in Table II. The percentages of each class at the various 
five levels are given in the second part of the table. 

The fluctuation introduced at high values of P brings in sources only at the 
places entered as F in Table II. The precise contribution of the fluctuation is 
uncertain within a factor 2. For curve II of Fig. 4, sources of class I would 
increase the total count at log P=1-5, and at log P=1-75 by about 30 per cent, 
while the count at log P=2-0 would be approximately doubled. 

Optical photographic magnitudes, uncorrected for spectrum red-shift effects, 
and making no allowance for galactic obscuration, are attached in brackets to 
the entries in the second part of Table II for the columns log P= —0°5, 0, 0:5, 
I, I1°5, 2. Since optical and radio magnitudes vary with distance in the same 
way when spectrum red-shift effects are omitted (these are in any case small for 
the radio magnitudes) the optical magnitudes change by 1:25 when the log P 
changes by 0-5. And assuming that all radio sources are associated with galaxies 
of the same intrinsic optical luminosity, the optical magnitudes vary as 2-5 log L 
for the various radio classes. Since log Z changes by o-5 from one radio class 
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to the next, the optical magnitude therefore also changes by 1-25 as we pass from 
one radio class to the next at fixed P. 

Remembering that Table I refers to class II sources (when the scale of P is 
10° w.m~*(c/s)-!), we note that the red-shift for log P=1 is close to o-15 
for aclass II source. At this red-shift, Humason, Mayall and Sandage give about 
+185 for the photographic magnitude of the brightest galaxy of the cluster 
0025 +2223. This value has been taken as the standard in Table II, other values 
being determined by differencing in the manner described in the previous 
paragraph. 

The magnitudes in Table II are therefore chosen as bright as seems reasonably 
possible. ‘They assume that all radio sources are associated with galaxies of 
abnormally great optical emissivity. While it is true that known optical 
identifications support this assumption, it is obvious that a selection effect exists 
in favour of making identifications with galaxies of abnormal luminosity. Also 
we may note that cluster 0025 + 2223 is at comparatively high galactic latitude 
b=—40°, where galactic absorption is only about o-1 mag. In attempting 
identification with radio sources it will usually be necessary to work at lower 
galactic latitudes where absorption is greater. Hence it does not seem likely 
that we have set the optical magnitudes too faint in Table II. 

The situation concerning the possibility of making optical identifications 
stands out clearly in the table. A satisfactory identification demands a highly 
accurate radio position, say +0°5’ in both coordinates, and so far it has not 
been found possible to obtain such positions except for particularly bright sources. 
And the remarkable point emerges from Table II that the bright radio sources 
are just the ones that contain the highest proportion of classes 1V, V and VI, 
all of which are distant objects associated with faint optical magnitudes. 

A thoroughgoing attempt to secure optical identifications has recently been 
made by Bolton and Minkowski, who have considered about 100 bright sources 
with flux levels above about log P= 1-5 (160 Mc/s). The attempted identifications 
with galaxies brighter than + 19-5 were successful in about 45 per cent of cases, 
which is about the percentage indicated by Table II, if the contribution F of 
class I is omitted. Table II suggests that some 35 per cent will be likely to fall 
between +17 and +19:5, as compared with the 30 per cent actually identified. 
The contribution F might have been expected also to yield some 30 per cent, 
whereas about 15 per cent were found at optical magnitudes bright enough to be 
associated with the nearby class I objects. 

The values of Table II do not hold out much hope that many further optical 
identifications will be forthcoming, unless either it proves possible to work to 
magnitude +20 (or even +21 if spectrum red-shift effects become important) 
or accurate radio positions become available down to log P~1. 

It is our impression that the following argument has surrounded the whole 
question of optical identifications. Taking the bright radio sources, roughly 
half of them cannot be identified with galaxies brighter than +19. Accordingly 
it seems fairly certain that these radio sources are distant objects. ‘Then how 
much more distant must the faint radio sources be? Our table shows that the 
last step of this argument may well be quite erroneous. The effect of the red-shift 
(not in the spectrum now, but on the apparent bolometric magnitude) can produce 
a situation in which classes IV, V, VI dominate the count at high, and even 
medium, values of P, but in which the intrinsically faint classes I and II entirely 
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dominate the count at low flux values. This effect has already been noted by 
Mills (4). 

It is of interest to consider the relation of our luminosity function to the 
recent work of Allen, Palmer and Rowson, on the surface temperatures of radio 
sources*. For the purpose of relating our results to the observations, we note 
that the L value for class VI is close to that of Cygnus A. Write T,,, for the 
measured surface temperature of Cygnus A Then if all radio sources possess 
the same physical size we expect the temperatures of our six classes to be given by 


(I, II, If, IV, V, VI) = (3.10%, 10~%, 3. 10-%, 10-7, 3. 10", 1)T¢,,. (27) 


Some uncertainty is attached to the appropriate value of T,,,. In the 
general surface temperature experiments, in which Cygnus A was treated in the 
same way as other sources, a value 7,,,~8x10°°K was found. Special 
information from the fuller analysis of Jennison et al. suggests a lower value of 
Toyg= 3 X 10° °K, however. A satisfactory compromise appears to be to accept 
the lower value of T,,,, but to continue to ignore the effect of the red-shift on 
the surface temperature, as we have already done in (27). The red-shift lowers 
the surface temperature by (1+2)-*. If we restrict ourselves to sources that 
make the important contributions to Table II, the values of z are never greater 
than 0-3 to 0-4, so that the (1+2)~* factor is not larger than the measure of 
uncertainty in T,,,. For the Jodrell Bank experiments, log P is on the average 
near 1-2, and at this value of P the red-shift correction is small for sources of 
classes I, II and III, although it amounts to a factor ~ 37-1 for class V. 

We set the surface temperatures in relation to Cygnus A as follows 


~ ~ ~ 7 ~ 7 « ~ 8 ~ 
T, = 10%, Ty, 23 x 10%, Ty,;210’, Tiyv=3 x10’, Tyz10°, Ty, 23 x 10°. 


(28) 


If the sources are not of the same physical size, if the intrinsically weaker sources 
are smaller in their dimensions, then the values of T for the lower classes have 
been set too low. ‘Turning to the percentages given in Table II for log P=1-25, 
we see that the most frequent source is of class III with surface temperature 
~ 10’ °K and that classes II and IV with temperatures ~3 x 10% and 3 x 10° 
are symmetrically distributed about ~10’°K. ‘This result is in close accord 
with the work of Allen, Palmer and Rowson, as also is the 20 per cent of sources 
at temperatures of order 10° °K. 

These considerations, both of the optical identifications and of the surface 
temperatures of radio sources, show that we have been realistic in the absolute 
values of L that have been associated with the various classes, and hence that 
our taking 10-*° w.m~? (c/s)~! as the unit of P is not an improper choice. 

As we have already mentioned, the mixing of the various luminosity classes 
results in a reduced slope for the log NP??—log P curve. The values given in 
the first part of Table II are plotted in Fig. 5. The slope of the log N—log P 
curve is now no greater than —1-6 in the important range of P, although there 
is a steepening to —1°7 at log P~1°5. The effect of the fluctuation at log P> 1-5 
has not been shown in Fig. 5. This would lift the curve in the region 1-5 < log P< 2, 
probably yielding a sub-Euclidean slope in this range, exactly as in the case of 
curve II of Fig. 4. 


* We are indebted to the Jodrell Bank workers for valuable discussions on their observations. 
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The present results are of course related to our choice, K(q)«q*, for the 
dependence of the probability function on age over the relevant generations 
(5th to gth or roth). A steeper rise of K would lead to a steeper log N—log P 
curve. For this reason no particular significance we feel can be attached to an 
amount of the order of o-1 or 0-2 in the slope. At the outset we noticed that 
whereas the latest Cambridge results yielded — 1-8, Mills has given reasons for 
believing that the slope should not exceed —1-5. In these circumstances our 
value of — 1-6 does not seem unsatisfactory. 

Although we have lost steepness in the slope for 0-4<logP<1 we have 
gained in the behaviour of the log NP®?—logP curve for logP<o-o. This 
section is now very close to the ultimate slope suggested by Ryle to be the most 
probable interpretation of an analysis by Hewish (of deflections rather than 
source counts) using the method of Scheuer. 

9. The count from individual lattice units—We return to a consideration of 
the luminosity class II alone. The specification of a red-shift z defines a distance 
cH-1z. Within this distance, if z is sufficiently large, there will be a particular 
generation m for which the volume of a unit cube, /,3, is comparable with 
4n/3(cH—z)%. Observation to red-shift z is expected to include the eight lattice 
members (or part of them) belonging to such a unit cube. The question arises 
as to how many sources are to be associated with each of these eight members 
and how does this number vary with z, and hence with P. An answer to this 
question can be given to within an uncertainty of a factor of about 2. ‘The 
uncertainty arises from connecting ideas that relate to the discrete model with 
values given by the continuous theory. ‘To remove the uncertainty it would be 
necessary to work entirely in terms of an explicit discrete model. (Such a project 
is feasible if automatic computing were used. If the present theory turns out to 
survive any further immediate attacks we shall consider it worth while to abandon 
altogether the simple methods of the continuous theory and to go over entirely 


to a discrete situation. ) 
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Referring back to (23), the number of sources counted to red-shift z,,, is 
given by integrating the function 


ray Kwesp (~3x) (29) 


over the appropriate portion of the (y,z) plane shown in Fig. 5. It will be 
recalled that yH~—! is the age of a source at the moment of emission or radiation, 
viz. x=7H—In(1+2z) where 7 is the present age. The relevant area of the 
(x,2) plane is that lying below the curve 

kz 
a (1+2)’ 
and enclosed by the z-axis, and by z=z,,. 

It is easy to estimate what fraction of the total count to red-shift z,, is 
contributed by the ‘‘last generation’’, this being defined as the last range 
of 4 in the interval of y. The fraction is obtained by integrating the same function 
(29) over the portion of the (y,z) plane enclosed by y=In(z,,/1+<2,,)—4, by 


x=In (30) 


“ml 


2=3,,, and by the curve (30). The fractional contribution is thus the ratio of 
the two integrals of (29) over the areas shown in Fig. 5. Introducing g=exp x 
again, we require the ratio of the following integrals 


z 


- 2m 22 dz 24 i 2) K(q) . 
| F aa) kz 7 dq, (31) 











kz 
2m 22dz 4'(1+2) K(q) 
paraclongie dq. 2 
; coil "i st 





a 
1+z k 

Computations for 0-2<2z<0-4 for K(q)«qg' and for k/z,'=15, give ratics 
close to 2/9. 

We turn now to the actual numbers N given in Table I. We have seen that 
if these numbers are interpreted per steradian they are in reasonable relation to 
the actual counts. Thus ~47N(z) is the order of the total count on the whole 
sky to red-shift z. We now see that ~2/g9 of this must be attributed to the 
“last generation’’, viz. ~ 3.N(z). 

The next step brings in the uncertainty of connection to the discrete model. 
We associate the ‘‘last generation’’, as defined above for the continuous 
picture, with the situation in the discrete model in which the observer sees the 
eight lattice units of generation m. Then each lattice unit of the “‘last 
generation ’’, contributes a total of ~ 3N(z) sources. In view of the uncertainty 
of association we can simply say that the lattice units making the maximum 
contribution to the count to red-shift z contain of the order of N(z) sources, 
where the functional dependence of N on z is given in Table I. (The present 
discussion essentially determines the order of the absolute value of the probability 
function K.) 

Reference to Table II shows that we have made important use of class II 
sources down to about logP=—o-5. ‘Table I shows that such a value of P 
corresponds to z~ 0-6, and at this red-shift N(z)~2500 sources. It follows that 
the maximum population of the oldest lattice units (those of the ~ roth generation) 
is some 2500 sources of class II. The maximum populations required for the 
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other luminosity classes are easily estimated from Table II. For example, the 
relative proportions of classes I and II at logP=-—o-5 is about 3/1. Our 
considerations down to this flux level therefore require the presence of some 
7500 sources of class I. Results for the various classes are as follows: 


Maximum number of sources required per lattice unit 


Class I 7500 
Class II 2500 
Class III 250 
Class IV 25 
Class V IO 
Class VI <1 


Since the number of galaxies per lattice unit is ~10° it is clear that the 
maximum populations in the oldest lattice units only become at all comparable 
with the number of galaxies for the cases of classes I and II. Only these cases 
require further consideration. 

Class II sources are comparable to the source Hydra A, for which Burbidge 
(13) has calculated a minimum total energy requirement of ~10°° ergs. At an 
emission rate of 4x 10"%ergsec~! (the value arrived at in Section 6), such an 
energy reservoir would last for about 2-5 x 10% sec. We require the main radio 
source activity throughout a lattice unit to last for an interval of time of the 
order of }H~-!. Thus if each of the ~ 10° galaxies were to become a radio source 
on one occasion during the relevant generation the number of sources existing 
at any particular moment would be of order 2:5 x 10!°/(4H~) x 10°.. H~! must 
be expressed in seconds, viz. H~-!'~3 x10!’ sec. The number is thus very 
close to 2500. 

A similar calculation for class I sources yields an equally satisfactory answer. 
Since the emission rate for class I is lower by a factor ~ 3 the time for which a 
given energy reservoir can maintain the radiation rate is increased by 3, and the 
number of sources existing at any moment of time is expressed by the same factor. 

It is of course the case that we do not literally require every galaxy to become 
a radio source. It would be equally sufficient if a sub-class of galaxies became 
radio sources on an appropriate number of occasions, say 10 per cent of galaxies 
ten times more frequently. Energy requirements do not forbid such a situation. 
Thus the gravitational energy of a galaxy of mass M and radius r is of order 
GM?/r; and for M=3x10"%0, r=15 kpc, this is ~5 x 10° erg. Hence .an 
energy reservoir ~ 10°° erg could be provided on a number of occasions. 

It is entirely possible that the classes [II-VI are simply cases where galaxies 
of abnormally large mass are involved. With M increased by ~ 10, gravitational 
energies are increased by ~ 102, which is of the order of difference between II 
and VI. The lower proportions required for the very strong radio sources 
could be simply a consequence of the small fraction of galaxies that possess masses 
of order 3 x 10/0. 

10. The distribution of sources on the sky.—The discrete model not only 
introduces fluctuations in the number counts, of the sort discussed in Section 7, 
but it also introduces fluctuations from strict isotropy in the distribution of 
sources on the sky. The discrete lattice units introduce patchiness. We have 
now to examine whether this patchy characteristic is consistent with the 
observations or not. 
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We note that MSH have drawn attention to a possible departure from isotropy, 
and that more recently Mills, Slee and Hill (14) report that between declinations 
—20° and — 50° they find a marked deficiency of sources in R.A. 00" to R.A. 022, 
(We have taken Mills’ sources and have divided them into 12 groups. Each group 
covers 2" in R.A. and the whole range of declination. We find only a negligible 
probability that the grouping is random. A test applied to the postulate that the 
source density varies in the ratio 4/3 from one area of sky of order 10* square 
degrees to another yielded a reasonable result.) 

We proceed now to the theoretical discussion, and we do this, first, for a 
single luminosity class. When a limiting flux P is given, the generations of our 
model can be ordered in the following way. There will be a particular generation, 
n(P) say, of which the observer will expect to see one complete unit cube with 
its associated 8 lattice units. For the n(P)+s generation (s=1, 0, —1, —2,...) 
the observer will expect to see ~8exp(-—s) lattice units. The series of values 
of s must be taken only for s<1, since for s>2 the formula derived from the 
continuous theory ceases to have any relevance. 

In the simple cubic model the 8 lattice units for s=o present the smallest 
total solid angle to an observer at equal distances from each of them. With / as 
the lattice spacing, the observer’s distance from the vertices is (+/3/2)/. Since 
the radii of the units (in our model) is 4/, the solid angle presented by each unit 
is 0-485 ster. The 8 units together therefore cover a fraction 0-31 of the whole sky. 

For the generation s = — 1, each lattice unit is smaller by e~"*. Thus a lattice 
unit of generation s= —1 at the same distance / subtends a solid angle smaller 
by ~e-?8. There are, however, e times more units from generation s= —1 
than from generation s=o. Generation s= —1 therefore subtends a total solid 
angle that is greater by at least e"’®. (Not all the units of s= —1 are as far away 
as /, although the majority of them are at distances of this order.) Hence 
generation s= —1 covers at least a fraction 0-31e"%=0-43 of the whole sky. 
In a similar way, one can easily see that generation s = — 2 covers at least a fraction 
0:3 1e?* =0-60 of the whole sky. 

The situation is somewhat more complicated for generation s=1. A lattice 
unit of this generation, if completely within range of observation, must subtend 
so large a solid angle that a simple inverse square law cannot be used. The most 
distant sources of the 8 lattice units of s=o, for an observer at the centre of a unit 
cube, demand observation not merely to 4/3//2 but to /(1/3+4+/3/2). If the 
sources of an s=1 dattice all lie within this distance, then the centre of the 
lattice unit cannot be further away from the observer than 


(; + 3) — ies 
3 2 3 
the radius of the unit being 4/exp (4). Hence the minimum solid angle subtended 
by such a unit must be 1-43 ster, if the unit is to be completely observed. Three such 
units would therefore cover at least a fraction 0-342 of the whole sky. 

It is of course the case that the formula 8exp(—s), derived from the 
continuous theory, does not require the value ~ 3, given by s=1, to be made up 
of three complete lattice units of generation s=1. The three could be made up 
by partial contributions from more than three lattice units. Indeed a system of 
partial contributions is certainly much more probable. Such a system would 
have a smaller departure from isotropy, since there would be more than three 
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patches on the sky, and each patch would contribute fewer sources. We shall 
keep, however, to the worst case in which the observer sees three complete lattice 
units of generation s=1. Summarizing we have: 


Solid angle subtended 


Generation Number of by each lattice unit Minimum fraction of 
s lattice units (ster.) the sky covered 
I 3 1°43 0°34 
° 8 0°48 o°31 
-I 22 0°25 0°43 
—2 60 O13 0°60 


It is useful to consider departures from isotropy by dividing these four 
generations into two pairs s=—1, —2; s=1,0. The first of these pairs has so 
many lattice units that it is essentially free from ambiguities; and the sum of 
the fractions of the sky covered adds almost exactly to unity, which is a convenient 
simplification. Since successive generations are expected to be anticorrelated we 
consider the patches of s= —1 to interlace those of s= —2. 

Write X for the total number of sources contributed by s= —1, —2, —3,.... 
The work of the previous section suggested that ~2/g X are contributed by 
$=-—1, about 2.2 by s=—2, etc. So many lattice units are contained in 
$= —3, —4, that we regard their contribution as being isotropic. ‘That is to say, 
49X/81 sources are distributed isotropically. It is easily seen that the average 
source density per unit area for a patch from s= —1 to that for a patch from 
s=—2 18 


[e604 + 5 (o43) | / [ 604+ xt (057) | =§, (33) 


where, for simplicity, we have slightly reduced the fraction of the sky covered 
by the s= —2 patches from 0-60 to 0°57. 

To test the effect of such a patchiness we took 27 squares arranged in three 
rows of 9. The squares were of two categories A and B. Equal probabilities 
were attached to each individual square being of type A or of type B. Decisions 
were made using random number tables: it turned out that 11 squares were 
assigned to type A and 16 totype B. The next step was to divide each such square 
into 100 sub-squares. A probability $p was attached to there being a source in 
each sub-square belonging to type A, and p was attached to there being a source 
in each sub-square belonging to B. The value of p was chosen so that a total 
number of sources closely equal to the 3C catalogue was obtained, and the 
probabilities were again controlled by random number tables. As a final detail, 
when a source happened to fall in a sub-square, its position in that sub-square 
was also randomized. The result is shown in the second section of Fig. 6, the 
first part being a plot of the 3C sources. 

Since the ratio 4/3 used in this test represented a greater measure of anisotropy 
than was estimated above in (33), it seems unlikely that the anisotropy arising 
from generations s=—1, —2 can bring us into a conflict with observation. 
It is of course true that statistical effects would be weakened if the total count X 
were increased—the ‘‘signal’’ carried by the probability ratio would stand out 
more in relation to the statistical ‘‘noise’’ as X increased, the signal-to-noise 
ratio improving essentially as X12. According to the recent Cambridge results, 
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about 8000 sources would be counted on the whole sky in a survey down to 
P~2x10-%w.m~*(c/s)"*. Of these, about 4000 would belong to the pair 
s=0, 1 which has still to be considered. The appropriate value of X for the 
pair s= —1, —2 would thus be about 4000 in such a survey, as compared with 
about 1000 used in Fig. 6. Hence an improvement in “signal-to-noise ratio ”’ 
by ~2 would be expected in such a survey, if all sources belonged to a single 
luminosity class. ‘The mixing of sources of different luminosity classes reduces 
the signal-to-noise ratio, however, first by reducing the value of X appropriate 
for each class, and second by increasing the ‘‘noise’’. The ‘‘noise’’ is 
increased because the different luminosity classes give different values of n(P). 
The distribution of patches is accordingly different for the different classes 
(the flux level P being specified). 


90°, 
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We turn to the pair s=1,0. Write Y for the total number of sources from 
the whole series s=1,0, —1,.... Then about 2/9 Y come from s=1, about 
14/81 Y from s=o, and 49/81 Y from s=1, —2,... (i.e. XXO°6 Y). Roughly 
half of the sources therefore arise in the pair s=1, o. 

Since the number of patches from both s=1 and s=o together is only of 
order 10, it is very reasonable to regard the positions on the sky of the patches 
as being anticorrelated to those of the s=o patches. If the two generations 
together are regarded as covering the whole sky then no crucial issue arises. The 
two generations cover approximately equal fractions of the sky, while they 
contribute numbers of sources in the ratio 9/7. 
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This ratio may somewhat underestimate the anisotropy, since a corresponding 
calculation using the discrete model gives a rather greater anisotropy. ‘The above 
ratio 9/7 is derived from the value 2/9 Y for s=1, 14/81 Y for s=o, obtained in 
Section g from the continuous picture. Instead, we can argue, perhaps more 
simply, that each lattice unit of s=1 contributes e‘* times more sources than 
each lattice unit of s=o. With three of the former and eight of the latter, the 
ratio is 3e“/3/8=1-43/1. But even with a ratio as great as 3/2 the fluctuation of 
source area density would merge with the smaller anisotropy of the generations 
s=—1, —2. A combination of the pairs s=—1, —2; s=1,0 could hardly 
yield a patch ratio significantly worse than 4/3 even if the two pairs were not 
anticorrelated. 

The serious question is whether the generations s=1, 0 really do cover the 
sky, and it appears to turn very largely on this issue as to whether a detectable 
anisotropy can be expected to arise or not. According to the above estimates, 
generations s=1, 0 cover 2/3 of the sky, but this estimate was very definitely a 
minimum value. The following argument is of interest in showing how far the 
fraction 0-31 for s=o is a minimum. 

It will be recalled that the value 0-31 was obtained for an observer at equal 
distances from all eight lattice units of generation s=o. A change of the observer 
to other positions increases the total solid angle subtended by these units. 
A calculation based upon a cubic lattice, of the effect of placing the observer in 
other positions, would attach too great an importance to the precise structure of 
the cubic lattice however. As stated in Section 4, one must take care not to use 
the cubic lattice model in such a way as to obtain results that are peculiar to this 
particular lattice, since in actuality the lattice may be expected to have no particular 
regularity. We therefore consider the following argument from which all lattice 
regularity is removed. 

Write v for the number of condensation centres per unit volume in generation 
s=o. 'To obtain eight such centres we require a sphere of radius a where 


4ma®v wil 


3 


Suppose a lattice unit at distance a subtends a solid angle 2, and approximate 
by taking the solid angle at distance r (<a) as Q,(a/r)*—once again this can 
only underestimate the solid angle. Then, on the average, the total solid angle 
subtended by the eight lattice units inside the sphere of radius a is 


a 
grvOet | dr = 47vQ, ~ 24 Q, 
0 


which is three times the value given by placing the whole eight units at equal 
distances a from the observer. 

This argument must not be given too great weight, however, since it considers 
no correlation at all in the positioning of the centres of the lattice units, but 
it does show how far varying the observer’s position, and varying the regularity 
of the lattice structure, could increase the total solid angle. To lift our 0-31 of 
Table IX to ~0-5 we need a much smaller increase than the factor 3 given by this 
argument. 
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Turning now to s=1, we need to increase the 0-34 of Table [X also to ~o°s. 
This could be achieved in two ways: by viewing more than three lattice units 
in a partial way—the advantages of this from the point of view of diminishing 
the anisotropy have already been mentioned—or by slightly decreasing the 
distance of the observer from the centres of each lattice unit. The estimate 0-34 
was based on a maximum distance consistent with complete visibility of all the 
sources of the three lattice units of s=1. ‘The solid angle is very sensitive to the 
precise distance of the observer, a slight decrease of distance producing a large 
increase of solid angle. ‘To increase the fraction of the sky covered by the three 
lattice units from 0°34 to ~0°5 it is only necessary to decrease the distance of the 
observer from their centres by some 20 per cent. 

We may sum up the above discussion as follows: dividing the generations 
contributing to the counts in the manner described, the generations s= — 3, —4,... 
are sensibly isotropic. It is unlikely that any anisotropy in the pair s= —1, —2 
can readily be detected. And if the pair s=o,1 fill the whole sky, it is also 
unlikely that anisotropy in s=o, 1 can be at all easily detected. If, however, 
there is a failure to fill the sky, a fluctuation of the source density over large 
areas of the sky might well be expected. ‘The relevant areas are of order 1 steradian. 
In this connection it is of interest that the fluctuation reported by Mills is over 
an area of this order. It is possible that Mills’ result arises from some systematic 
error, but if so the error would illustrate the experimental difficulty of giving a 
decisive answer to whether such large-scale fluctuations exist or not. 


APPENDIX 

It might appear at first sight that results of the previous section can 
be reproduced by introducing an age—luminosity function L(7) instead of the 
age—probability function K(7). As mentioned earlier, the two cases are not 
equivalent. While A(7) affects the integrand in (19) the function L(r) affecis 
the limits. It is proposed to give an example where we take K(r)=const. and 
introduce L(7). We will consider the log N —log P curve on the hypothesis that 
L(r) «x exp (3H7) on analogy with the form chosen for K(7). 

The count to power level P is given by 


N(P)= | lr exp (—3Hx) dx dz (34) 


where the area of integration is specified in Fig. 7 by the shaded region. The 
upper curve « is given by 


kz 
—] — = 
x n (1+2)2,’ (35) 
and the lower one f by 
x =41n [P2?(1+2)7], (36) 


where x is measured in units of H-'. 

Let (2, x2), (%3, x3) denote the intersections of « and 8; 2,>23. If P is 
large enough, x,;<0, and we have to integrate from 2; to z,. If P is small, so 
that y;>0, we have to integrate from 2,/k—2, to 24, and from 2, to 22, where 2, 
is the point where the curve f intersects OZ, as shown in Fig. 7. In the absence 
of a power-age correlation, the curve 8 would simply be a straight line parallel 
to the x-axis. 
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It is best to express N, P as functions of the parameter 2, and to obtain their 
values for different values of z.. From the figure it can be seen that P has a 
maximum value P,, when the curves «, 8 touch one another. ‘This happens 
when z,=%3;=0°25. For P>P,, N=o. Here again we encounter the effect of 
discreteness on the N —P curve. 

The double integral can be evaluated for various values of z,. Some values 


are given in the following table, with z,'/k=1/15 as before. 
Ze Zs Bq N P NP?!? 
0°25 0°25 —_— ° ‘O81 ° 
O'5 orl! —_ 0°26 0658 4°38 
I'O — “086 2°618 03125 14°46 
1°s — 123 67526 ‘O154 12°47 
2°0 _— 163 15°203 “0082 II‘29 
% 


(2X5 ) 







(x) 
x=Ln ae (6) 


(1+z)z, ~\ 





x= 3 Ln [P2%() +2] 








Fic. 7 


Here also NP®? increases initially as P decreases. But the increase is very 
rapid and is maintained only over a small range in P (a range ~2), while the 
decrease is gradual. It therefore does not reproduce the observed results. An 
exponential Lo exp (4H7r) would give similar behaviour, the case Loc exp (37) 
was adopted because of the resulting simplicity of the integrals involved. 

Thus though the amount of radiation contributed by a cluster of age 7 is 
the same when L(r)=f(7), K(r)=const. as when K(r)=f(7), L(7)=const. 
the two cases are not equivalent as far as the N— P curve is concerned. 

This discussion, however, does not rule out the possibility of a combination 
of a power-age correlation with an age-dependent probability. 


St John’s College, Fitzwilliam House, 
Cambridge: Cambridge : 
1961 Fune. 
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EXTRAPOLATION OF THE NUMBER-FLUX DENSITY RELATION 
OF RADIO STARS BY SCHEUER’S STATISTICAL METHOD 


A. Hewish 
(Received 1961 June 14) 


Summary 


An attempt is made to assess the usefulness of Scheuer’s statistical method 
of analysing the records derived from a phase-switching receiver. A Monte 
Carlo technique is described which enables the probability distribution of the 
recorded deflections to be computed rapidly in EDSAC for any assumed 
distribution of radio sources. Some typical models are considered and it is 
concluded that, subject to a small enough receiver noise, the statistical method 
gives information about the source distribution at flux densities considerably 
lower than those at which the sources may be counted individually. The 
method is applied to some observations at 178 Mc/s using the method of 
aperture synthesis. It is shown that the log N—log S relation cannot have a 
uniform slope of —1°5. A slope of —1°8 yields good agreement provided 
that the source density is suitably truncated at small flux densities. 





1. Introduction 

Records obtained during a survey of radio stars may be analysed in two 
different ways. The usual method is to extract the positions and flux densities 
of individual sources; in this case the limiting flux density reached in the survey 
is determined either by the sensitivity of the receiver or by confusion errors 
arising from the presence of more than one source in the reception pattern 
simultaneously. To reduce confusion errors to an acceptable level, several 
authors have suggested that it is necessary to limit the number of sources 
extracted to about one source in 20 beam areas. 

A different analysis of the data may also be carried out in which no attempt 
is made to isolate individual sources, but which gives statistical information 
concerning their number and flux density. In this method the record is sampled 
at equal intervals of time without regard to the occurrence of particular sources. 
For the case of an interferometer the deflection D at a given instant is defined 
as the modulus of the quasi-periodic record, while for a total power system it 
is simply the total deflection as shown in Fig. 1. When a sufficient number of 
samples has been obtained, the result may be presented as the distribution 
function P(D) where P(D) dD is the probability of a deflection lying in the 
range D to D+dD. ‘The shape of P(D) depends, of course, on whether an 
interferometer or a total power system is used and in this paper only interferometric 
systems are considered. 

The quantity D has the dimensions of flux density S, and the physical 
meaning of P(D) is clarified by considering the ideal case of an aerial which 
receives with uniform sensitivity over a prescribed beam area. In this case, 
for sufficiently large values of D corresponding to sources of high flux density 
which are clearly resolved, we have P(D)=P(S) where P(S) is the chance that 
a source of flux density S to S+dS lies within one beam area; for smaller 
values of D this equivalence no longer holds, since there is then a finite chance 
of several sources lying in the beam area at once, and P(D)->o as D-o while 
P(S) may increase without limit. When the reception pattern has an arbitrary 
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(b) 


Fic. 1.—Values of D derived by sampling the output of (a) a phase-switching receiver and 
(b) a total power receiver. 


profile P(D) and P(S) are similarly related but a constant of proportionality 
is involved. 

The precise relation between P(D) and P(S) has been fully discussed by 
Scheuer (1957). He showed that it was not possible to derive P(S) uniquely 
from P(D): what may be done, however, is to assume some relation for P(S) 
which enables P(D) to be calculated theoretically and then compared with the 
observed curve. If the experimental P(D) curve agrees with the theoretical 
one, then the assumed P(S) is consistent with the data. A population of radio 
sources is usually defined by N(.S) where N(S) is the total number of sources 
per steradian having flux densities greater than S; but 


N(S)oc | * P(S) dS 
S 


so that assumed models of N(S) may be directly related to P(S). By this 
process it is possible, in principle, to extrapolate the observed N(S) relation 
towards low flux densities where the sources can no longer be counted with 
accuracy owing to confusion errors. ‘The utility of the statistical method in 
any practical application will, however, depend on the sensitivity of P(D) to 
variations in the assumed values of N(S). 

In this paper an attempt is made to assess the usefulness of the statistical 
method. Analytical methods which allow P(D) to be computed when (5S) 
is given have been described by Scheuer (1957). The computations are 
laborious, however, and do not lend themselves readily to rapid comparison of 
different N(.S) models. A method has therefore been adopted in which 
P(D) is computed by a Monte Carlo technique using the electronic computer 
EDSAC II. The basis of the method is to consider a typical beam area and to 
populate it, with appropriate statistical fluctuations, according to some assumed 
N(S) model. ‘The different contributions to the total deflection are then 
combined with random phase and the computation is repeated until P(D) 
may be determined with sufficiently small statistical uncertainty. Using 
pseudo-random numbers generated in the computer it is possible to derive 
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P(D) for a given model in a matter of minutes. The computational method 
is given in Section 2, and in Section 3 the general utility of the statistical method 
is discussed in the practical case where receiver noise must be taken into account. 

Finally, in Section 4, the method is applied to results obtained at 178 Mc/s 
in a survey using the principle of aperture synthesis. The radio observations 
have been fully described in another paper (Scott, Ryle and Hewish 1961) 
where it was explained how the observed P(D) relation was automatically 
derived from the computer as part of the routine synthesis programme. 
A summary of the available data concerning the observed N(S) relation 
obtained during this and other surveys has also been given (Scott and Ryle 
1961) and it has been shown that effects due to the angular diameter and 
clustering of sources can only be relevant to a small fraction of the total number 
(Leslie 1961). In the present paper the statistical method is used to extrapolate 
the observed N(S) relation to smaller flux densities. It is shown that the shape 
of the observed P(D) relation cannot be explained by any N(S) model such that 
N(S)oc. S~} for all values of S. Models with N(.S)oc S~!* are in good agreement 
with the observations, provided that the source density is appropriately 
truncated at small flux densities. 


2. Computation of P(D) by a Monte Carlo method 


For a phase-switching interferometer the value of D at any given sampling 
instant is the resultant of a large number of contributions of different amplitude 
and phase arising from the sources lying within the beam area at that moment. 
The essence of the Monte Carlo method is to consider a typical beam area, to 
populate it randomly with sources whose mean distribution is defined by some 
N(S) model and to sum the different contributions with random phase. The 
computation is repeated for a large number of hypothetical beam areas, the 
population being varied randomly from sample to sample, until the theoretical 
P(D) may be derived with a reasonable statistical certainty. About 3000 samples 
were usually taken which gave statistical errors comparable to those of the 
observations. The computation will now be discussed in detail. 

(i) The correction for the experimental reception pattern.—The basic compu- 
tation was carried out for an assumed ideal reception pattern having uniform 
sensitivity over the entire beam area. With a realistic reception pattern a source 
of given flux density may give rise to a variety of deflections according to 
its position in the beam. ‘This effect could have been introduced into the 
computation without difficulty by means of another parameter subject to random 
variations, but it involved less computing time to adopt an alternative procedure 
in which the N(.S) model was suitably modified as a preliminary step. Since 
the sources lying within each interval of flux density S to S+dS will give 
deflections of different magnitudes appropriate to a random sampling of the 
reception pattern, the result obtained when using a graduated reception pattern 
is exactly simulated by smoothing the N(.S) model. ‘The smoothing process 
is an exact convolution if the reception pattern is formulated as the probability 
of obtaining a given deflection when a constant source is randomly positioned 
within it. ‘The essential randomness is then introduced in calculating the 
random population of the beam area as discussed below. When calculating in 
this way it is assumed that both the phase and amplitude of the quasi-periodic 
record may be considered as independent random variables. This assumption 
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is justified if, as in the present case, the interference pattern contains many 
cycles. 

The equivalent areas of the experimental 178 Mc/s reception pattern, 
considered constant over small intervals of relative sensitivity, are shown in 
Fig. 2; these values were derived from the published reception pattern (Scott, 
Ryle and Hewish 1961) and include the effect of side lobes. The logarithmic 
intervals adopted are the same as those employed in the main computation 
discussed below. 


Beam 
area, 
ino 


200 


100 











Relative sensitivity 


Fic. 2.—The experimental 178 Mc/s reception pattern expressed as increments of beam area for 
equal logarithmic intervals of relative sensitivity. 


The result of convolving N(.S) with the probability distribution of the 
reception pattern may be denoted by N(D) where N(D) is the total number 
of deflections > D which combine to give the net deflection appropriate to each 
beam area. It should be noticed that whereas N(S) is expressed as the number 
of sources per steradian, N(D) is the number of deflections per beam area. 
When N(S) obeys a simple power law, e.g. N(S)oc.S~'5, N(S) and N(D) are 
related by a simple scale factor; but when N(S) is terminated, or suffers a 
change of slope, the corresponding change in N(D) will be smoothed by the 
convolution process. An example of this is shown in Fig. 3 where N(D) is 
given for an abruptly terminated N(S) model. 

(ii) Populating the sample beam area.—For simplicity the full range of D 
was partitioned into discrete intervals AD such that log,, AD/D=1/10. Within 
a given partition the AN deflections were all assumed to have the same magnitude. 
The first stage of the main computation was to populate the partitions randomly 
according to an average distribution defined by some particular model N(D). 
This process is equivalent to constructing a typical beam area containing a random 
distribution of sources of different flux density. ‘The manner in which the 
required population was computed depended upon the value of AN and 
four different methods were employed. 

(a) AN<1. This case corresponds to the small chance of obtaining a large 
deflection in one beam area. Here it was sufficient to generate a random number 
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in the range o-1000 and to add a deflection if the random number fell inside 
a given range of 1000xAWN integers. Different selections of integers were 
chosen to represent the different partitions and one trial was allowed for each 
partition. This method will lead to errors, however, if there is a finite chance 
of obtaining more than one deflection in a given partition during a run of, 
say, 3000 sample beam areas, that is, when 3000 AN?~1 or AN~1°8x 107%. 
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Fic. 3.—N(D) derived by convolving the N(.S) model (c) depicted in Fig. 5 (b) with the reception 
pattern as expressed in Fig. 2. 


(b) 10 ?<AN<s5. ‘To achieve the correct population when a given partition 
may contain zero, or more than one deflection, a range of partitions was 
considered such that SAN~10. In a true random population AN will vary, 
from sample to sample, by an amount 6N such that 5N?~ZAN where 6N is 
approximately normally distributed. For each sample beam area a deviation 
5N was generated by taking a random number and using it to choose one term 
in a table of random deviates. This determined the value TAN’=LAN+6N 
for the sample and a further XAN’ random numbers were then generated 
and used to distribute the LAN’ deflections amongst the different partitions 
with a probability proportional to AN. 

(c) 5<AN<s50. For small enough values of D there is a negligible chance 
of obtaining zero deflections within any partition. In this case it was sufficient 
to compute a series of values AN’=AN+8N in order to specify the random 
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population. One random number was therefore generated for each partition 
and used to assign a random deviation as in case (b). 

(d) AN>50. When the number of deflections in any partition is sufficiently 
large, corresponding to many weak sources in the beam area simultaneously, 
there is no longer any necessity to compute the random walk summation by 
Monte Carlo methods. The contribution from any one partition will be a vector 
of random phase whose amplitude is governed, from sample to sample, by the 
well-known Rayleigh distribution function. For many of the N(D) models a 
simple power law was adopted in the region of small D and in this case the 
contribution due to the Rayleigh component was easily calculated. If 


N(D)=AD", where A is constant, the mean square deflection d(D®) due to 
components in the range D to D+dD is given by 

d(D®)=nAD"* dD 
so that 


D?= | * nAD"™ dD 
Dy 


nN \"? 
Deus = (- as) Dy. 


giving 





To compute the precise value of the deflection in each sample beam area 
a random number was generated and used to assign a deviation drawn at random 
from a table of deviations appropriate to a Rayleigh distribution having a 
standard deviation given by Dpys/ 1/2: 

For some of the models N tended to a constant value for small values of D 
and Dpys had to be derived by numerical integration. 

(iii) The random walk summation.—Having populated a sample beam area 
it remains to add the component deflections with random phase. For each 
component a random phase angle was generated by choosing a unit vector from 
a table of vectors tabulated at 3° intervals in the range o-27, the choice being 
determined by a random number. Throughout the computation random 
numbers were generated as required by the multiplicative process M,,,,=M.M,, 
(modulo 2°) with M=51!5 and M,=1. Pseudo-random series of this nature 
are non-repetitive within a period of 2°? (Moshman 1954) and an initial test 
was carried out to check the randomness in a run of the required length. As 
an additional safeguard an arrangement was made whereby the calculation was 
halted if M,=M; no such event occurred. 

The resultant of the random walk summation for each sample beam area 
was accumulated as a histogram, together with the total population of each 
partition. 

(iv) Receiver notse.—In any practical determination of P(D), deflections 
arising from noise in the receiving system will be present in addition to those 
due to signals from the aerials. In the method of aperture synthesis the output 
of the receiver is recorded on punched tape, and the final integration, normally 
performed by the output filter of a receiver, is carried out by a convolution 
process in the computer itself. The equivalent filter is arranged to have a narrow 
frequency response centred on the mean frequency of the quasi-periodic signal 
and this ensures that noise fluctuations alone will produce a quasi-periodic 
output whose modulus is specified by a Rayleigh distribution. The deflection 
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due to receiver noise has a character which is thus exactly equivalent to case 
(i(d)) and it may be treated in precisely the same manner. Putting 
Days =(D? + Dy*)¥*, where D,? is the mean square deflection due to noise 
alone, makes allowance for both quantities simultaneously. 

(v) The validity of the Monte Carlo method.—Because a number of 
approximations have been used in the computation it is desirable to test the 
validity of the method by using it to reproduce well-known results before 
proceeding to unknown cases. ‘T'wo such tests were applied. In the first a 
model was taken such that N(D)ocD~!> and the random population was 
assigned as described above. In the random walk summation, however, the 
individual deflections in all the partitions were given equal magnitudes. In this 
case P(D) is the well-known Rayleigh distribution and it is seen from Fig. 4 (a) 
that the Monte Carlo result is in very good agreement. 


P(D) 





eT ee ee ee ee ee ee a ee 


D 





Fic. 4 (a).—The continuous curve shows a Rayleigh distribution. The histogram indicates the 


same result obtained by the Monte Carlo method. 











Fic. 4 (b).—P(D) derived analytically by Scheuer for N(S)& S-*5 (continuous curve) compared with 
the histogram given by the Monte Carlo method. 
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As a second test the same model was taken and the computation carried out 
in the normal way. For this model, corresponding to a uniform spatial 
distribution of radio sources, P(D) has been computed analytically by Scheuer 
(1957) and both results are shown in Fig. 4(b); the agreement is again seen 
to be close. It should be noted that “curve fitting’’ plays no part in the 
Monte Carlo method. A given model of N(S), and hence of N(D), defines 
the scale of D absolutely. 


3. The value of the statistical method 


To estimate the usefulness of the statistical method it is necessary to investigate 
the sensitivity of P(D) to changes in N(D). Only when such changes are greater 
than the statistical uncertainty inherent in a practical determination of P(D) can 
the method be usefully applied. Initially it is assumed that the receiver is ideal, 
so that noise components do not contribute to P(D). Two types of model have 
been considered, as shown in Figs. 5 (a) and 5 (b), such that N(S)ocS~1° and 
N(S)oc.S~''8; the latter model is more relevant to the actual distribution of 
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Fic. 5 (a).—Models in which N(S)S7*. Fic. 5 (b).—Models in which N(S)XS-}. 


sources. In both cases the effect of arbitrarily terminating N(S) at various 
values of S was investigated. P(D) corresponding to the different cases is plotted 
in Figs. 6(a) and 6(6); for slope—1-5 a significant change in P(D) can be 
detected when N(S) is terminated at S=o-o1, corresponding to the level at 
which there are 50 sources per beam area. For’slope—1-8 there is a more 
pronounced change in P(D) when N(S) is terminated at S=o-o1, where there 
are 250 sources per beam area. 

It is therefore clear that, in the absence of receiver noise, the statistical 
method gives information about N(S) at flux densities considerably lower 
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than those at which the sources may be counted individually. To investigate 
the effect of receiver noise the model shown in Fig. 5 (6) was taken and P(D) 
computed for two different RMS noise levels. The results are shown in 
Figs. 7(a) and 7(b). While P(D) is modified, as would be expected for 
sufficiently large noise levels, it is seen that the sensitivity of P(D) to changes 
in N(S) is not appreciably diminished provided that the maximum of P(D), 
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Fic. 6 (b).—P(D) derived from the N(S) models illustrated in Fig. 5 (b). 


in the absence of receiver noise, occurs at a flux density somewhat higher than 
the RMS noise level. This result is to be expected in view of the Rayleigh-like 
character of P(D). When the noise level is high enough to cause a significant 
displacement of the maximum of P(D), but not so high that the shape of P(D) 
is dominated by it, the receiver noise must necessarily be accurately known if 
the statistical method is to be applied without loss of precision. 

The discussion of the preceding sections has been confined to interferometric 
systems. For a total power recording the situation is similar but the maximum 
of P(D) is then located at a flux density near the mean d.c. level of the receiver 
output. While it is simple to compute P(D) for a given N(S) model by the 
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Fic. 7 (a).—P(D) derived from the N(S) models shown in Fig. 5 (b) with an assumed RMS receiver 
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Fic. 7 (b).—The corresponding P(D) curves with a receiver noise of 0 
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method described in Section 2, except that the contributions to D are added 
directly and not with random phases, it is difficult in practice to define a flux 
density corresponding to the mean receiver output. The comparison of theoretical 
and observational P(D) relations must then rely on the shape of the curves and 
the position of the maximum cannot be used as a criterion of agreement. 


4. Application of the method to observations at 178 Mc/s 
As described in detail in another paper (Scott, Ryle and Hewish 1961) 
the reduction of the observations made at 178 Mc/s using the method of aperture 


synthesis included the automatic computation of P(D). 
experimental curve derived from a survey of 0-3 steradians 


In Fig. 8 is shown the 
(a = 20h 40m cade 1g! rS™, 


5=40°— 44°); the RMS noise level was determined in a separate experiment and 


found to be 0-085 x 10°76 w.m~? (c/s)"1. 
from this survey, and from “‘ whole sky”’ surveys usin 
« bd . 


Counts of individual sources derived 


different arrangements 


of the 178 Mc/s interferometer, gave the relation N(.S) oc S~ 18 for sources having 


flux densities greater than 2 x 102° w.m~? (c/s)~1 (Ryle a 
In attempts to account for the observed P(D) curve i 


nd Clarke 1961). 
n terms of appropriate 


models, attention was concentrated in the region of the maximum of P(D), 


and models were sought which fitted this portion of the curve most closely. 
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For larger values of D the statistical uncertainty is relatively greater since P(D) 
is small; moreover, we have 


[ P(D) dDxN(S) 


e 


(cf. Section 1) for large D so that no information is gained by studying P(D) in 
this region, where N(S) is already known from the observed source counts. 


P(D 
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Fic. 8.—P(D) derived from observations at 178 Mc/s. 


As a first attempt to account for the observations, models were taken such 
that N(S)ocS-1%, The three models are illustrated in Fig. g(a) and the 
corresponding P(D) curves are given in Fig. g(6). It is immediately clear 
that none of the models is satisfactory. From the trend of the curves it can 
be seen that the maxima are either too high (model (c)) or occur at too large 
a flux density (models (6b) and (a)). Even without evidence provided by 
individual source counts it may be concluded that no model having N(S) oc. S~15 
can give a satisfactory explanation of the observed N(S) relation. 

In a second attempt to find a suitable model the observed N(S) relation was 
extrapolated at constant slope (— 1-8) and successively truncated at a series of 
flux densities. These models were discussed in a different context in the previous 
section and are shown in Fig. 5 (6). The corresponding P(D) curves are 
shown in Fig. 7 (a) and it was found that model (e) fitted the observed P(D) 
within the statistical error. For this model the truncation occurs at a flux 
density of 0-35 x 10°76 w.m~? (c/s)~! corresponding to the level at which one 
source occurs,,on the average, in two beam areas. 

It was also possible to account for the observations by extrapolating from 
a flux density near the limit of the observed number counts with a slope of 
—1-5 and truncating at a flux density of 0-18 x 10°-°® w.m~* (c/s) "1. Models 
of this type and the corresponding P(D) curves are indicated in Figs. 10 (a) 
and ro (5b). 
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Fic. 9 (a).—Three models such that N(S)&S-!* used in an attempt to explain the observed P(D). 
N(S) obtained directly from source counts at 178 Mc/s is denoted by the thick line. 
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Fic. 9 (5).—P(D) computed from the models shown in Fig. 9 (a). The observed P(D) is 
also shown (thick line). 
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Fic. 10 (a).—N(S) models in which the observed source counts are extrapolated with a slope 


of —1°5. 
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Fic. 10.(b).—P(D) computed for the models shown in Fig. 10 (a). 
(thick line). 
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The observed P(D) is shown 
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Fic. 11 (a).—Models in which the observed source counts are extrapolated with a slope of —1'9. 
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Fic. 11 (6).—P(D) computed for the models shown in Fig. 11 (a). 
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Since a sudden truncation corresponds to an unrealistic physical situation 
another series of models was taken in which the \(S) relation was extrapolated 
to zero flux density with a slope of —1-0. This series is shown in Figs. 11 (a) 
and 11 (d) and it was again possible to derive a curve in reasonable agreement 
with the observations. 
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. Fic. 12.—Extrapolations of the observed source counts which are consistent with the 
observed P(D) relation, 

The three models which fitted the observations most closely are shown 
together in Fig. 12 and it appears that any satisfactory model must lie within 
the boundaries indicated. It is interesting to notice that the number 
of sources begins to depart appreciably from the expectations of a linear 
extrapolation of the observed source count at a flux density in the range ovr 
to 1:0x 1077 w.m-*(c/s)"1.. In view of the strong evidence favouring an 
extragalactic origin of the radio sources (Ryle and Clarke 1961) this decrease 

i is presumably a cosmological effect. The discontinuities in the curves will, of 
i course, be smoothed by any spread in absolute luminosity. The magnitude of 
this effect may be estimated with reference to the smoothing action of the 
reception pattern already discussed in Section 2(i). Fig. 3 indicates the 
rounding-off which would occur for a sudden cut-off in the source distribution ; 
if the spread in absolute luminosity were of the order of 100: 1. 
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THREE SOUTHERN SPECTROSCOPIC BINARIES 
William Buscombe and Pamela M. Morris 
(Received 1961 June 19) 


Summary 


Radial velocities from 51 spectra, obtained in 1955-60, are presented 
for three single-lined binaries, HD 184035, 114911 (7 Mus) and 217792 
a PsA). HD 184035 is a newly discovered binary. For the other stars early 
observations at the Lick Chile Station and at the Cape have been included 
in the determination of the orbital elements, which for all three stars are 
listed in Table VI. 





Introduction.—A programme is under way at Mt Stromlo to determine the 
radial velocities of certain bright stars in the fundamental catalogue N30 (1), 
for which reliable values were not available. In addition to the intrinsic interest 
in solving the orbits of binary stars, there is the hope that ephemerides of their 
velocity variation may be useful for standardizing measures from future objective- 
prism surveys. 

Most of the observations on this programme were made with the 3-prism 
spectrograph at the Cassegrain focus of the 30-inch reflector. For a few of the 
spectrograms we are especially indebted to other observers, who have observed 
the stars at our request during particular phases, as indicated in ‘lables II] and V 
below. A few plates were also taken with the 2-prism spectrograph at the New- 
tonian focus of the 74-inch reflector. For convenience, details of the various 
instruments from which observations have been included in this discussion are 
listedin l'ableI. ‘The three stars are named in Table II, of which the final column 
indicates the spectrographs from which observations are available. 


TABLE I 
Instruments 
Observatory Telescope Spectrograph dispersion at Hy 
Royal Cape, South Africa 24-inch refractor 4-prism 19 A/mm 
Lick, Santiago, Chile 36-inch Cassegrain 2-prism 20 
I-prism 37 
Mt Stromlo, Australia 30-inch Cassegrain 3-prism 36 
74-inch Newtonian 2-prism go 
Tas_e II 
The binary stars 
HD Star Mag. Spectrum R.A. (1950) Dec. Obs. 
114911 7 Mus 4°79 B8 V 13%rr™-8 = =—67°37’ L1, St2, St3 
184035 201G Sgr 5°90 Aj Ill Ig 30 ‘7 —4009 St2, St3 
217792 a PsA 5°13 Fo IV 23 00 °7 —35 OI C4, L2, St3 
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In an earlier communication to the Society (2), the authors have reported on 
the progress of the general programme, with special reference to experience in 
measuring the spectra of standard-velocity stars and in the selection of wavelengths 
for the absorption lines. Further details are now in the press, including some 
results of the extension of this work to the Newtonian spectrograph. 

Since the number of observations of each star is not large, or uniformly 
distributed through the different phases, the method outlined by Irwin (3) was 
followed. After the period is found, points on the assumed velocity curve are read 
at certain fractions of the amplitude K. Using these values in various combina- 
tions, auxiliary diagrams are drawn to determine the most probable values of e 
and w by the use of Irwin’s tables. Under the circumstances, the computation of 
adjustments to the preliminary orbital parameters seems unwarranted. 

HD 114911.—From the very approximate value of its radial velocity as deter- 
mined at the former Lick Southern Station (4), some authors have considered 
H{D 114911 as a possible member of the Scorpio-Centaurus Association. How- 
ever, the luminosity class indicates that it is a foreground object. ‘The extrema 
in velocity observed in 1924 and in 1960 establish the orbital period, but the 
number of measures secured may be regarded as the fewest possible for computing 
a reliable orbit. 


TaBLe III 
Observations of HD 114911 

JD Phase V O-C Notes 
re km/sec 
6979'221 O°412 — 3°5 +3°3 Underexposed 
6982°228 562 + 30°7 =a Newtonian spectrogram 
7002°179 *560 + 30°3 +0°4 Overexposed 
7064°995 "700 +44°8 —o'2 
7076992 "299 —44°0 —4°4 
7082'936 ‘597 +28°6 —7°7 Weak comparison 
7098'903 *3905 —37°3 = 
7108'934 ‘896 + o'2 —4°7 Observer: A. Przybylski 
7II1'g29 ‘046 —53°7 +1°7 Observer: H. Gollnow 
7111°937 ‘046 —45'8 +98 Underexposed 
7130°871 0°993 — 39°5 “#9 


The velocities measured at Mt Stromlo are listed in Table III, with the values 
of the residuals from the spectrographic orbit. ‘The computed velocity curve is 
shown in Fig. 1. The parameters derived are shown in Table VI. ‘The absorp- 
tion lines are rather broad, and Huang (5) has estimated the rotational velocity 
component in the line of sight as 240 km/sec. 

HD 184035.—The spectrum of HD 184035 shows unusually sharp lines for an 
A-type star, from which the velocity measures are remarkably consistent. As 
exposures of more than one hour became necessary with the Cassegrain spectro- 
graph on the photographic emulsion 103aO available in 1960, the last three 
observations were obtained with the 74-inch reflector and the Newtonian spectro- 
graph. 

Since few cycles of the velocity variation have been observed, the period cannot 
be determined with high accuracy. However, Fig. 2 shows that the observations 
listed in ‘Table IV are well represented by the velocity curve computed from the 
orbital parameters given in Table VI. 
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Fic. 1.—The radial velocity curve of y Mus. Open circles represent Lick observations, 
filled circles Mt Stromlo. 


TaBLe IV 
Observations of HD 184035 
Jb Phase Vo O-C Notes 
7 km/sec 
6739°165 0°942 +56°5 +15°1 
6778-065 an +21°1 + 5°5 
7065°284 "454 —22°4 ~~ 29 
7077°245 "040 +72°3 — 86 
7099°197 *787 —30°8 + 12 
7115147 "235 +57°0 — 6:0 
7122°138 ‘747 —44°7 - 09 
7172°988 "742 —43°0 + 2°1 
7174°964 ‘169 +86°1 + 4°0 Underexposed 
7177°947 ‘810 —24°2 - oF 
7182°940 893 +14°7 — 13 
7210'915 "942 +43°6 + 2°2 
7216°899 236 + 68:0 + 5°! 
7237°917 “780 —31'8 + 2°4 Newtonian spectrogram 
7242°905 859 — 36 — 21 Newtonian spectrogram 
7245°892 0°505 — 36:2 — 02 Newtonian spectrogram 


HD 217792.—Before we made our observations it was informally suggested 
to the authors that this star might be a cepheid. In view of the long period, 
the low luminosity and the shape of the velocity curve, this seems most unlikely. 
It is noteworthy that several of the earlier measures, especially from the Cape 
Observatory (6), deviate systematically from the adopted curve for our recent 
observations. If a secular change in the orbital elements has occurred, larger 
values of both e and w would satisfy the early data somewhat better than the 
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Fic. 3.—The radial velocity curve of 7 PsA. Open circles, Lick; 
filled circles, Mt Stromlo. 
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TABLE V 


Observations 

Vo 
km/sec 
+ 84 
— 92 
+ 27°5 
+29°6 
+28°2 
— 78 
—14°1 
—14°1 
+ OF 
—I10°! 
—10°6 
—12°0 
—12°0 

fohze) 
— 05 
+ 1§°S 
- 71 
— 6°9 
= 259 
Pr E°% 
+20°7 
+15°4 
+ 12°0 
—11°8 


of HD 217792 


O-C Notes 


Observer: G. Hagemann 
—I'l Observer: 5. Gaposchkin 


—o'l Underexposed 


+2°9 Observer: A. Przybylski 
+3°1 Observer: H. Gollnow 


+19 Observer: H. Gollnow 


TaBLe VI 


Orbital parameters 


114911 
20°0052 
—S81 
56°5 

o'12 

120 
20606°722 
15°4 
0°36 


154035 217792 
4°625 356°567 
+119 —9°9 
72°5 21°5 

0°09 0°70 

315° 10 
36702°433 34968-2066 
4°6 75°3 

ors O13 


The early observations from the Mills Expedition of Lick 
Observatory (4), especially those near the sharp velocity maximum, were useful 


The Stromlo velocities are listed in ‘Table V, and 
* all the available data are shown in Fig. 3. A few spectrograms of this star have 


1 been re-measured, for confirmation, by Mr John Graham and by Miss Willy van 
\ Pelt. The orbital elements are listed in ‘l'able VI. 


ise : Mount Stromlo Observatory, 


Canberra, A.C.T., 


Australia: 


1961 Fune. 
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Summary of Edinburgh paper 


SPECTROPHOTOMETRIC MEASUREMENTS OF EARLY-TYPE 


STARS* 


VII. ReEsuLts AND DiscussION FOR 10 Stars OF MK Type Bs5 
AND 7 STARS OF MISCELLANEOUS ‘TYPES 


H. E. Butler and G. I. Thompson 


Summary 


This paper concludes the spectrophotometric programme on early-type 
stars. Itisin two parts: 


SECTION I 

Equivalent widths are given for 56 absorption lines in the spectra of 10 
stars of MK type B5. The data are strictly comparable with those of previous 
papers of this series. Comparing their behaviour with that of the B3 stars, 
no unexpected trends have shown up except that He 1 6678 and 5675 have 
become appreciably luminosity sensitive. Al 111 has disappeared while S 11 
and Fe 11 are perceptible. 

There is a suggestion of an emission band in the supergiants centred on 
4745A. This is tentatively identified as C11. There is also a broad but 
unidentified absorption in 67 Oph at 4272-80A. 

Unusually sharp emission wings are visible in the profiles of Ha and H/ 
in the supergiant 5 Per. An appreciable change in the absorption component 
of Ha has also been observed. 


SECTION 2 


This section deals with stars that did not fit into the previous parts, or 
whose spectra have changed. 

One Ao supergiant and one Ao dwarf have been reduced in a manner 
comparable with the rest of the programme. 

Three Be stars, previously described and showing Ha in emission, were 
later re-observed and showed Ha in a changed form. New reductions show 
that the changes were, in fact, confined to the hydrogen lines and to the 
D-lines of sodium. ‘This latter observation suggests that the D lines are 
largely circumstellar rather than interstellar. 

Further plates of w Ori (B3 IIIe) confirm the existence of an unexplained 
depression centred at 4348 in the red wing of Hy. 

The unusual and nearly featureless spectrum of the star 59 Cyg is studied 
and its relation to abnormal stars such as ¢ Per and ys Per commented on. 





